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Abstract 


Ultra- Wideband  (UWB)  impulse  radar  is  inherently  a  noise-limited  technology.  While  sub¬ 
nanosecond  pulses  achieve  good  range  resolution,  the  resultant  extreme  bandwidth  exacts  a 
toll  on  system  performance  in  terms  of  noise  power  entering  the  receiver.  Data  must  be 
processed  efficiently  to  yield  an  acceptable  signal  to  noise  ratio  (SNR),  and  to  enhance 
structural  details  of  the  return  pulse.  This  is  achieved  by  a  combination  of  ensemble  averaging 
to  reduce  thermal  fluctuations  and  the  newly  developed  Background  Noise  Conditioning 
(BNC)  method,  a  statistical  signal  processing  technique  used  to  remove  isolated  spectral 
interferers.  BNC,  developed  at  DRDC  Ottawa,  periodically  examines  spectral  noise  properties 
to  automatically  configure  an  appropriate  notch  filter. 

For  the  radar  data  obtained  in  this  investigation,  averaging  over  a  relatively  small  number  of 
pulses  is  insufficient  to  raise  the  desired  returns  appreciably  above  the  noise  floor.  Typical 
SNR  values  for  one  received  pulse  are  around  18  to  19  dB.  Beyond  ten  pulses,  SNR  improves 
approximately  linearly.  Integration  of  about  80  pulses  results  in  a  SNR  maximum  of  about  35 
dB ,  with  an  associated  improvement  of  ~ 15  dB  over  the  untreated  data  (that  is,  for  situations 
where  BNC  is  not  applied).  For  pulses  >  80 ,  the  averaging  takes  place  over  a  long  enough 
time  period  that  returns  become  blurred  due  to  target  motion.  This  effectively  defines  the 
upper  limit  of  pulse  returns  that  can  be  averaged  without  additional  processing  to  focus  the 
returns.  Part  of  this  report  is  devoted  to  modeling  and  characterizing  the  impulse  shape  and  its 
associated  spectrum.  As  well,  the  minimum  detectable  target  velocity  is  extracted  from  data 
and  compared  with  theory,  showing  good  agreement. 

Wall  effects  have  been  investigated — there  was  an  experimental  determination  of 
approximately  8dB  attenuation  losses  due  to  a  cindercrete  wall.  The  same  wall,  once  soaked 
with  water,  caused  only  an  additional  2dB  loss  (both  of  these  are  one-way  measurements).  A 
model  for  wall-induced  clutter  is  developed,  providing  expressions  for  the  pulse  return  shape 
and  duration. 

The  imaging  of  objects  buried  in  snow,  using  both  ID  ranging  and  2D  SAR  techniques 
establishes  that  objects  can  be  located,  although  the  current  resolution  is  such  that  it  would 
make  identification  challenging.  A  major  concern  is  the  processing  technique  itself,  which  is  a 
form  of  area  Moving  Target  Indicator  (MTI),  where  returns  from  stationary  targets  are 
partially  removed  in  the  processing  stage.  These  results  will  be  used  to  develop  algorithms  for 
extracting  stationary  targets  without  apriori  knowledge  of  the  background. 
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Resume 


De  par  sa  nature,  le  radar  a  impulsions  a  bande  ultra-large  (UWB)  correspond  a  une 
technologic  limitee  par  le  bruit.  Bien  que  les  impulsions  de  duree  inferieure  a  une 
nanoseconde  donnent  une  bonne  resolution  en  distance,  la  largeur  de  bande  extreme  qui  en 
resulte  nuit  aux  performances  du  systeme  en  termes  de  puissance  de  bruit  entrant  dans  le 
recepteur.  Les  donnees  doivent  etre  traitees  efficacement  pour  produire  un  rapport  signal/bruit 
(S/B)  acceptable  et  pour  rehausser  les  details  structuraux  de  l’echo.  A  cette  fin,  on  utilise  une 
combinaison  de  l’etablissement  de  la  moyenne  d’ensemble  afin  de  reduire  les  fluctuations 
thermiques  et  de  la  methode  recemment  elaboree  de  conditionnement  de  bruit  de  fond  (BNC), 
technique  de  traitement  statistique  des  signaux  utilisee  pour  eliminer  les  signaux  brouilleurs 
spectraux  isoles.  Le  BNC,  elabore  par  RJDDC  Ottawa,  examine  periodiquement  les  proprietes 
du  bruit  spectral  en  vue  de  configurer  automatiquement  un  filtre  coupe-bande  approprie. 

Pour  les  donnees  radar  obtenues  au  cours  de  cette  etude,  l’etablissement  d’une  moyenne  pour 
un  nombre  d’ impulsions  relativement  petit  ne  suffit  pas  pour  porter  les  echos  voulus 
notablement  au-dessus  du  bruit  de  fond.  Les  valeurs  types  du  rapport  S/B  d’une  seule 
impulsion  regue  sont  de  l’ordre  de  18  a  19  dB.  Au-dela  de  dix  impulsions,  le  rapport  S/B 
s’ameliore  de  fagon  a  peu  pres  lineaire.  L’integration  d’environ  80  impulsions  donne  lieu  a  un 
rapport  S/B  maximal  d’environ  35  dB ,  avec  une  amelioration  associee  de  -15  dB  des  donnees 
non  traitees  (c.-a-d.  pour  les  cas  ou  le  BNC  n’est  pas  applique).  Pour  plus  de  80  impulsions, 
l’etablissement  de  la  moyenne  a  lieu  pendant  une  periode  assez  longue  pour  que  les  echos 
soient  brouilles  a  cause  du  deplacement  de  la  cible.  Cela  definit  effectivement  la  limite 
superieure  des  echos  dont  un  peut  etablir  la  moyenne  sans  traitement  supplementaire  pour 
cibler  les  echos.  Une  partie  de  ce  rapport  est  consacree  a  la  modelisation  et  a  la  caracterisation 
de  la  forme  d’ impulsion  et  du  spectre  connexe.  De  plus,  la  vitesse  minimale  detectable  de  la 
cible  est  extraite  des  donnees  et  comparee  avec  la  valeur  calculee  :  il  y  a  une  bonne 
concordance. 

On  a  etudie  les  effets  de  paroi  —  une  experience  a  permis  de  determiner  un  affaiblissement 
d’environ  8  dB  du  a  une  paroi  en  beton  de  machefer.  La  meme  paroi,  trempee  d’eau,  causait 
un  affaiblissement  supplementaire  de  seulement  2  dB  (dans  les  eux  cas,  il  s’agissait  de 
mesures  unidirectionnelles).  Un  modele  est  elabore  pour  le  clutter  du  a  la  presence  d’une 
paroi,  permettant  d’exprimer  la  forme  et  la  duree  des  echos. 

L’imagerie  d’objets  enfouis  dans  la  neige,  au  moyen  de  techniques  de  telemetrie  ID  et  de 
techniques  RSO  2D,  a  permis  d’etablir  que  ces  objets  peuvent  etre  localises,  bien  que  la 
resolution  actuelle  soit  telle  que  1’ identification  poserait  un  defi.  Une  des  preoccupations 
principales  est  la  technique  de  traitement  elle-meme,  qui  fait  appel  a  une  forme  d’ indication 
de  cibles  mobiles  (MTI)  de  zone,  et  dans  laquelle  les  cibles  fixes  sont  partiellement  eliminees 
au  stade  de  traitement.  Ces  resultats  seront  utilises  pour  elaborer  des  algorithmes  pour  extraire 
les  cibles  fixes  sans  connaissances  a  priori  du  bruit  de  fond. 
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Executive  summary 


Defence  R&D  Canada  (DRDC)  has  an  active  research  program  in  many  aspects  of  Ultra- 
Wideband  (UWB)  radar  systems,  including  simulation,  data  analysis  and  hardware 
prototyping.  UWB  impulse  (or  short-pulse)  radars  are  capable  of  producing  high-resolution 
images  by  combining  the  performance  benefits  of  sub-nanosecond  impulse  technology  and 
synthetic  arrays.  One  particular  application  of  this  is  in  the  area  of  through- wall  sensing, 
where  relatively  low  center  frequencies  enable  good  radio  frequency  (RF)  penetration  of 
typical  building  materials.  Another  application  involves  the  use  of  UWB  radar  systems  for 
search  and  rescue  of  avalanche  victims.  This  is  a  Canadian  multi-department  venture 
including  Parks  Canada,  the  National  Search  and  Rescue  Secretariat,  and  the  Department  of 
National  Defence. 

This  report  outlines  the  current  target  discrimination  capabilities  of  UWB  impulse  radar 
systems  to  operate  through  simple  building  materials  and  snow.  While  sub-nanosecond  pulses 
achieve  good  range  resolution,  the  resultant  extreme  bandwidth  exacts  a  toll  on  system 
performance  in  terms  of  noise  power  entering  the  receiver.  Impulse  radar  lacks  the  filtering 
and  down-conversion  of  narrowband  systems,  and  individual  returns  are  severely  impaired  by 
thermal  noise,  clutter  and  external  interference.  Consequently,  UWB  impulse  radar  is 
inherently  a  noise-limited  technology. 

To  compensate,  data  must  be  processed  to  yield  an  acceptable  SNR,  and  to  enhance  structural 
details  of  the  return  pulse.  This  is  achieved  by  a  combination  of  ensemble  averaging  to  reduce 
thermal  fluctuations,  and  application  of  the  newly  developed  Background  Noise  Conditioning 
(BNC)  method,  a  statistical  signal  processing  technique  used  to  remove  isolated  spectral 
interferers.  BNC,  developed  at  DRDC  Ottawa,  periodically  examines  spectral  noise  properties 
to  automatically  configure  an  appropriate  notch  filter.  Another  attribute  considered  in  this 
report  is  the  measurement  of  target  velocity,  traditionally  obtained  in  narrowband  systems  via 
the  Doppler  frequency  shift  arising  from  target  motion.  Impulse  radar  systems  are  baseband, 
and  have  no  carrier  frequency  to  exploit.  However,  this  information  can  be  obtained  directly 
in  the  time-domain  by  measuring  changes  in  pulse  arrival  times. 

This  report  demonstrates  improved  target  discrimination  capability,  achieving  enhanced  signal 
to  noise  ratios  (SNR),  and  verification  of  target  velocity  thresholds  using  experimentally 
obtained  data.  These  results  will  be  used  to  develop  algorithms  for  extracting  both  moving 
and  stationary  targets  without  apriori  knowledge  of  the  background.  One  approach  could  be  to 
establish  a  data  bank  of  target  signatures  and  use  these  to  process  data  returns  via  correlation 
techniques. 
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Sommaire 


R  &  D  pour  la  defense  Canada  (RDDC)  a  un  programme  de  recherche  active  sur  bien  des 
aspects  des  systemes  radar  a  bande  ultra-large  (UWB),  y  compris  la  simulation,  l’analyse  des 
donnees  et  le  prototypage  materiel.  Les  radars  a  impulsions  (ou  a  impulsions  courtes)  UWB 
sont  capables  de  produire  des  images  a  haute  resolution  en  combinant  les  avantages  de 
performance  de  la  technologie  des  impulsions  de  duree  inferieure  a  une  nanoseconde  et  le 
radar  a  synthese  d’ouverture.  Une  application  particuliere  de  cette  technique  est  dans  le 
domaine  de  la  detection  a  travers  une  paroi,  ou  des  frequences  centrales  relativement  basses 
permettent  une  bonne  penetration  RF  de  materiaux  de  construction  types.  Une  autre 
application  concerne  l’utilisation  de  systemes  radar  UWB  pour  la  recherche  et  le  sauvetage  de 
victimes  d’ avalanche.  II  s’agit  d’une  entreprise  canadienne  pluriministerielle  comprenant 
Parcs  Canada,  le  Secretariat  national  de  recherche  et  sauvetage  et  le  ministere  de  la  Defense 
nationale. 

Le  present  rapport  decrit  les  capacites  actuelles  de  discrimination  de  cibles  des  systemes  radar 
a  impulsions  UWB  a  travers  de  simples  materiaux  de  construction  et  la  neige.  Bien  que  les 
impulsions  de  duree  inferieure  a  une  nanoseconde  donnent  une  bonne  resolution  en  distance, 
la  largeur  de  bande  extreme  qui  en  resulte  nuit  aux  performances  du  systeme  en  termes  de 
puissance  de  bruit  entrant  dans  le  recepteur.  Le  radar  a  impulsions  est  depourvu  du  filtrage  et 
de  la  conversion  vers  le  bas  des  systemes  a  bande  etroite,  et  les  echos  individuels  sont 
fortement  alteres  par  le  bruit  thermique,  le  clutter  et  le  brouillage  externe.  Par  consequent,  de 
par  sa  nature,  le  radar  a  impulsions  UWB  correspond  a  une  technologie  limitee  par  le  bruit. 

Pour  y  remedier,  les  donnees  doivent  etre  traitees  pour  produire  un  rapport  S/B  acceptable  et 
pour  rehausser  les  details  structuraux  de  l’echo.  A  cette  fin,  on  utilise  une  combinaison  de 
l’etablissement  de  la  moyenne  d’ensemble  afin  de  reduire  les  fluctuations  thermiques  et  de 
l’application  de  la  methode  recemment  elaboree  de  conditionnement  de  bruit  de  fond  (BNC), 
technique  de  traitement  statistique  des  signaux  utilisee  pour  eliminer  les  signaux  brouilleurs 
spectraux  isoles.  Le  BNC,  elabore  par  RDDC  Ottawa,  examine  periodiquement  les  proprietes 
du  bruit  spectral  en  vue  de  configurer  automatiquement  un  filtre  coupe-bande  approprie.  Un 
autre  aspect  considere  dans  le  present  rapport  est  la  mesure  de  la  vitesse  de  la  cible,  obtenue 
jusque-la  au  moyen  de  systemes  a  bande  etroite  grace  au  decalage  Doppler  produit  par  le 
deplacement  de  la  cible.  Les  systemes  radar  a  impulsions  fonctionnent  en  bande  de  base  et  ne 
peuvent  pas  mettre  en  valeur  une  frequence  porteuse.  Toutefois,  l’information  de  vitesse  peut 
etre  obtenue  directement  dans  le  domaine  temps  par  mesure  de  la  variation  du  temps  d’arrivee 
des  echos. 

Le  present  rapport  demontre  la  capacite  amelioree  de  discrimination  de  cibles,  l’obtention  de 
meilleurs  rapports  signal/bruit  (S/B)  et  la  verification  de  seuils  de  vitesse  des  cibles  au  moyen 
de  donnees  experimentales.  Ces  resultats  seront  utilises  pour  elaborer  des  algorithmes  pour 
extraire  les  cibles  fixes  sans  connaissances  a  priori  du  bruit  de  fond.  Une  des  approches 
possibles  consisterait  a  etablir  une  banque  de  donnees  de  signatures  de  cibles  et  a  utiliser  ces 
donnees  pour  traiter  les  echos  au  moyen  de  techniques  de  mise  en  correlation. 
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1.  Introduction 


Defence  R&D  Canada  (DRDC)  has  an  active  research  program  in  many  aspects  of  Ultra- 
Wideband  (UWB)  radar  systems,  including  simulation,  data  analysis  and  hardware 
prototyping.  One  particular  application  is  in  the  area  of  through-wall  sensing,  where  relatively 
low  center  frequencies  enable  good  radio  frequency  (RF)  penetration  of  typical  building 
materials.  A  second  application  involves  the  use  of  UWB  radar  systems  for  search  and  rescue 
of  avalanche  victims.  This  is  a  Canadian  multi-department  venture  including  Parks  Canada, 
the  National  Search  and  Rescue  Secretariat,  and  the  Department  of  National  Defence. 

UWB  radars  are  capable  of  producing  high-resolution  images  by  combining  the  benefits  of 
sub-nanosecond  impulse  technology  and  synthetic  arrays.  A  challenging  but  nonetheless 
necessary  requirement  for  military,  and  civilian  application  of  UWB  radar  is  the  ability  to 
detect  and  characterize  possible  targets  through  relatively  dense,  inhomogeneous  materials. 

Previous  work  in  UWB  impulse  radar  has  focused  on  modeling  and  simulation  [1-3]  and  on 
broad,  theoretical  system-level  aspects  [4-6],  rather  than  the  processing  of  experimentally 
obtained  data.  This  document  attempts  to  address  the  imbalance  by  undertaking  an 
investigation  of  how  best  to  process  and  extract  useful  information  from  experimental  data. 
The  original  impetus  for  this  investigation  was  to  examine  the  feasibility  of  standoff 
application  for  UWB  radar  surveillance.  Some  of  the  questions  addressed  are: 

•  How  easily  can  a  target  be  observed,  and  at  what  range? 

•  How  many  pulses  need  to  be  integrated  to  obtain  a  good  signal,  and  how  can  this  be 
quantified? 

•  Using  a  sphere  as  a  target  means  the  RCS  is  approximately  constant,  but  how  well  does 
the  Friis  transmission  formula  apply  for  impulse  signals? 

•  Given  the  variable  nature  of  the  signal  (i.e.,  non-sinusoidal),  how  can  we  best  determine 
performance  metrics  such  as  SNR? 

•  What  is  the  typical  noise  background  and  what  are  the  characteristics  of  the  RF  channel? 

•  What  is  the  spectral  response? 

This  report  outlines  progress  that  has  been  achieved  in  target  discrimination  capability, 
achieving  enhanced  signal  to  noise  ratios  (SNR),  and  verification  of  target  velocity  thresholds. 
After  a  brief  discussion  of  the  experimental  configuration,  the  following  subjects  are 
examined: 

Filtering  and  noise  reduction  —  UWB  impulse  radar  is  inherently  a  noise-limited 
technology.  While  sub-nanosecond  pulses  achieve  good  range  resolution,  the  resultant 
extreme  bandwidth  exacts  a  toll  on  system  performance  in  terms  of  noise  power  entering  the 
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receiver.  This  section  examines  ways  to  efficiently  process  the  raw  data  to  achieve  acceptable 
Signal  to  Noise  levels. 

Signal  characterization  —  Impulse  waveforms  are  baseband-there  is  no  dominant  frequency 
component  as  is  the  case  with  narrowband.  The  question  arises  then  on  how  to  best 
characterize  the  signal.  Since  the  very  act  of  transmission  distorts  the  signal  we  need  a  better 
understanding  of  its  properties  in  both  time  and  frequency  domains. 

System  performance  —  The  experimental  apparatus  described  in  this  report  is  not  as  much  a 
radar  system  as  a  collection  of  instruments  for  transmission  and  sampling  of  data. 
Nevertheless,  we  can  still  piece  together  some  of  the  operating  parameters  to  estimate  the 
performance  in  terms  of  signal  to  noise  ratio  and  propagation  phenomena. 

Target  velocity  processing  —  Impulse  radar  transmissions  are  carrier-free  and  therefore 
unable  to  exploit  Doppler  processing.  However,  target  motion  is  still  observable  by  examining 
changes  in  target  range  from  pulse  to  pulse.  The  basis  for  processing  now  becomes  the  change 
in  arrival  time  of  the  pulse.  In  effect,  we  employ  a  form  of  Area  Moving  Target  Indicator 
(Area  MT1 ),  comparing  the  complex  envelopes  of  successive  scans  to  identify  target  motion. 
This  section  develops  the  governing  equations  for  minimum  detectable  velocity,  and 
demonstrates  how  this  information  can  be  extracted  from  experimental  data. 

Wall  considerations  —  As  concrete  was  selected  to  be  the  building  material  of  interest  for 
this  project,  we  need  to  determine  the  effects  of  moisture  content  on  signal  propagation  and 
attenuation.  To  our  knowledge,  little  has  been  mentioned  in  the  literature  regarding  moisture 
content  and  to  what  extent  it  will  increase  attenuation.  Accordingly,  a  concrete  block 
(cindercrete)  wall  was  constructed  to  obtain  transmission  loss  data  for  both  dry  and  wet 
blocks. 

Additionally,  for  through- wall  imaging  the  building  material  acts  as  a  reflector,  with  its  own 
characteristic  radar  cross  section  (RCS).  This  section  describes  an  adaptation  of  altimeter 
radar  applied  to  wall-induced  clutter,  and  derives  expressions  for  the  shape  and  duration  of 
wall  reflections. 

Application  to  avalanche  victim  detection  —  In  any  two-  or  three-dimensional  radar  imaging 
system,  the  basic  element  of  target  information  is  a  time-delayed  impulse  response.  This 
section  first  examines  one-dimensional  range  data  to  measure  radar  returns  from  various 
objects  buried  in  a  mound  of  snow,  and  then  supplies  some  recent  results  from  2D  SAR 
processing. 
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2.  Data  collection 


A  series  of  experiments  were  carried  out  using  an  improvised  UWB  radar  system  (see  Figure 
1)  to  track  the  motion  of  a  metallic  sphere  as  it  rotates  about  a  fixed  point.  From  trial  to  trial, 
the  nominal  range  remained  fixed  while  the  motion  of  the  sphere  traced  out  a  circular  path  one 
meter  in  diameter.  Sets  of  data  were  recorded  for  the  following  ranges:  7,  14,  21,  28,  35,  42, 
and  49m.  Essential  to  any  imaging  application  is  the  collection  of  raw  data  for  further 
processing.  This  configuration  provides  a  means  of  generating  and  transmitting  short-time 
pulses,  collecting  the  return  from  a  reflecting  target,  and  associating  with  that  echo  a  return¬ 
time  delay.  The  resulting  data  was  collected  and  written  onto  a  disc  for  post-processing. 


AVTECH  HV1-C-T 1 


Target:  Rotating  Sphere 


Nominal  Range 


Figure  1.  Configuration  for  data  collection. 

The  data  was  collected  in  frames  of  2000  samples  each1,  with  sampling  increments  of  50ps 
(sampling  frequency  of  20GHz),  to  provide  a  frame  duration  of  100ns.  The  pulse  repetition 
frequency  was  set  at  1kHz2,  resulting  in  a  pulse-to-pulse  interval  of  1ms.  The  elapsed  time  at 
any  given  frame  is  therefore  k  x  1ms,  where  k  =  frame  number.  For  longer  ranges,  the 
oscilloscope  trigger  time  was  delayed  to  keep  the  returns  in  the  desired  range  gate.  A  typical 
target  return  is  shown  in  Figure  2.  Features  to  note  are  the  prominent  return  at  14m, 
representing  the  sphere,  and  another,  less  coherent,  but  strong,  signal  at  the  origin.  This  is 
leakage  from  the  transmit  antenna,  and  as  it  turns  out,  is  an  unavoidable  artifact  since  it  is  this 
transmitted  pulse  that  is  used  to  set  the  trigger  on  the  oscilloscope.  In  between  these  two 
returns  is  a  region  dominated  by  thermal  noise,  external  interference  and  clutter,  referred  to  on 
the  whole  as  background  noise.  Part  of  this  document  examines  ways  to  minimize  the 
influence  of  this. 


1  Tektronix  sampling  oscilloscope  model  TDS-7404. 

2  AVTECH  pulse  generator  model  AVH-HV1-C-T1 
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Amplitude  [V] 


Figure  2.  Target  return  at  14m. 


The  configuration  depicted  in  Figure  1  is  obviously  not  an  operational  SAR  system,  but  it  can 
be  considered  equivalent  to  a  single  element,  or  discrete  position,  within  the  synthetic  array. 
With  this  single-platform,  monostatic  arrangement,  we  are  effectively  reduced  to  a  one¬ 
dimensional  ranging  device,  since  there  is  no  way  to  distinguish  angles  of  arrival.  However, 
as  will  be  developed  later  in  the  report,  the  motion  of  the  target  can  be  used  to  our  advantage 
in  developing  an  Inverse  SAR  (ISAR)  processor. 
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3.  Filtering  and  noise  reduction 


UWB  impulse  radar  is  inherently  a  noise-limited  technology.  While  sub-nanosecond  pulses 
achieve  good  range  resolution,  the  resultant  extreme  bandwidth  exacts  a  toll  on  system 
performance  in  terms  of  noise  power  entering  the  receiver. 

As  is  typical  of  impulse  radar  returns  that  lack  the  filtering  commonplace  in  narrowband 
systems,  individual  returns  are  severely  impaired  by  thermal  noise,  clutter  and  external 
interference.  A  variety  of  noise  sources  impact  the  filtering  requirements  of  UWB  impulse 
radar  systems: 


•  out-of-band  interference 

•  in-band  interference 

•  thermal  noise 

Data  must  be  processed  efficiently  to  yield  an  acceptable  signal  to  noise  ratio  (SNR),  and  to 
enhance  structural  details  of  the  return  pulse.  This  is  achieved  by  a  combination  of  ensemble 
averaging  to  reduce  thermal  fluctuations  and  the  newly-developed  Background  Noise 
Conditioning  (BNC)  method,  a  statistical  signal  processing  technique  used  to  remove  isolated 
spectral  interferers  [7].  The  BNC  algorithm  periodically  examines  spectral  noise/interference 
characteristics,  and  automatically  configures  an  appropriate  notch  filter. 

3.1  Out-of-band  interference 

Interfering  sources  that  lie  outside  of  the  band  of  interest  can  be  removed  by  standard  filtering 
techniques  [8].  In  the  early  stages  of  this  investigation,  data  was  processed  offline.  In  that 
situation,  a  bandpass  filter  was  implemented  in  software.  Specifically,  a  15-tap  FIR  filter  with 
a  passband  of  1.5  to  3.0  GHz  was  applied.  The  frequency  and  phase  response  are  displayed  in 
Figure  3.  A  Hamming  window  was  applied  and  the  output  normalized  so  that  the  magnitude 
response  is  OdB  at  the  center  frequency. 
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More  recently,  hardware  was  purchased  from  Mini-Circuits  to  carry  out  band  pass  filtering  as 
data  is  collected.  A  practical  aspect  of  these  filters  (SHP-1000  and  SLP-2950)  is  the  SMA- 
type  connector,  allowing  easy  addition  or  removal  from  our  prototype  system.  The  measured 
response  is  shown  in  Figure  4a.  Two  sets  of  such  filters  were  combined  to  produce  a  cascaded 
system  as  displayed  in  Figure  4b. 
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Figure  4.  Mini-Circuits  measured  filter  response. 

a) ,  single  highpass/lowpass  combination; 

b) .  two  cascaded  highpass/lowpass  filters. 


3.2  In-band  interference 

Removing  in-band  interference  is  somewhat  more  problematic.  By  its  very  definition, 
wideband  transmission  allows  interfering  sources  to  enter  the  receive  chain.  Further,  the 
interference  exists  in  two  main  forms — relatively  narrowband  “discrete”  sources  and 
broadband  noise;  we  require  different  approaches  to  remove  these  impairments. 
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Consider  for  example,  the  frequency  spectrum  displayed  in  Figure  5.  This  data  was  obtained 
passively,  that  is,  the  radar  was  not  transmitting  at  the  time.  This  is  an  indication  of  typical 
noise  background3.  Interspersed  throughout  the  spectrum  is  a  range  of  relatively  large  spectral 
components,  each  of  which  contributes  to  the  overall  received  signal,  introducing  artifacts  in 
the  waveform. 


Figure  5.  Background  noise,  no  radar  signal  present 
(dc  component  removed,  fs  =  20GHz). 


3.3  Background  noise  conditioning 

To  better  gauge  the  effect  of  noise  and  interference  on  the  signal,  major  peaks  were  first 
identified  and  a  filter  was  designed  to  “notch-out”  the  undesired  components.  The  result  was 
the  removal  of  low-frequency  ripples  in  the  waveform,  leading  to  an  improvement  in  SNR. 

The  problem  with  this  approach  is  that  it  requires  human  intervention  to  identify  unwanted 
components.  As  we  have  no  reason  to  assume  background  noise  will  remain  stationary,  the 
notch  filter  constructed  at  one  moment  will  not  always  apply  the  next.  Even  with  a  signal 
present  as  in  Figure  6,  we  readily  observe  a  number  of  peaks,  but  most  are  out  of  the 
bandwidth  of  interest,  which  extends  up  to  approximately  4GHz.  Only  one  distinct  peak  at 
approximately  2.5GHz  is  visible.  Components  buried  within  the  signal  are  the  major  concern 
in  this  situation. 


3  Additional  samples  were  obtained  at  different  times  and  different  locations.  The  plot  shown  here  is 
qualitatively  similar  to  other  data  sets. 
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Figure  6.  Raw,  unfiltered  spectral  data  (dc  component  removed). 

To  address  this,  an  algorithm  was  developed  to  automatically  configure  an  appropriate  notch 
filter.  The  basic  principle  of  operation  is  that  relatively  few  spectral  components  will  exceed  a 
given  value,  and  will  be  observed  in  the  spectrum  as  spikes  (refer  again  to  Figure  5).  Since  we 
do  not  apriori  know  this  cut-off,  it  is  difficult  to  provide  a  specific  value  for  which  we  could 
consider  that  component  a  spike. 

A  systematic  approach  to  this  is  to  examine  the  statistical  properties  of  background  noise.  For 
example,  by  using  the  spectral  data  to  generate  a  histogram  showing  the  relative  frequency  vs. 
spectral  magnitude  (Figure  7),  one  can  obtain  an  estimate  for  the  probability  distribution  f(x j, 
of  the  noise.  Following  this,  the  method  mimics  standard  radar  detection  techniques  where  an 
amplitude  threshold  value,  £  distinguishes  noise  from  signal  power.  Empirically,  it  was  found 
that  acceptable  performance  resulted  when  the  cut-off  was  set  at  P(%)  =  0.97 ,  where  formally, 

<? 

p(g) = J  f(x ) dx ■  ( 1 ) 

o 

This  technique  is  called  Background  Noise  Conditioning.  The  variable  x  represents  spectral 
magnitude.  The  phrase  “acceptable  performance”  is  somewhat  subjective,  but  in  general  terms 
the  processed  signals  displayed  a  noticable  improvement  in  SNR  (>  6dB)  without  significant 
change  to  the  overall  structure  of  the  waveforms. 
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P  =  0.85 


spectral  magnitude 

Figure  7.  Typical  probability  density  function  f(x)  of  background  noise  frequency  spectrum. 


It  is  envisioned  that  in  operational  radar,  this  type  of  processing  could  be  carried  out  by 
periodically  suppressing  a  pulse,  with  the  receiver  listening  passively.  Over  time,  a  database 
could  be  established,  generating  a  set  of  best-fit  parameters  to  model  a  representative  pdf  such 
as,  say,  a  lognormal  or  a  Rayleigh  distribution. 


3.4  Thermal  noise 

Due  to  the  spectral  content  of  UWB  impulse  signals,  commonly  used  narrowband  filtering 
methods  are  not  helpful.  The  extreme  bandwidth,  together  with  a  lack  of  carrier  signal  means 
we  have  no  recourse  to  heterodyning  techniques  and  the  ensuing  intermediate  frequency  (IF) 
filtering.  To  obtain  a  useable  SNR,  the  system  performs  multiple  sample  integrations  to 
“average  out”  noise  fluctuations. 

This  technique  is  applied  to  UWB  processing  in  the  form  of  ensemble  averaging.  To  enhance 
signal  levels,  multiple  frames  are  collected  for  processing.  A  number  of  features  remain 
relatively  unchanged  across  frames.  This  could  include: 

•  Stationary  objects  in  the  illumination  area 

•  Effects  of  system  ringing  (for  example,  antennas,  cabling,  instrumentation  effects) 

•  Coupling  of  pulses  from  transmit  to  receive  antenna. 

Impairments  of  this  type  must  be  removed  from  the  raw  data.  In  a  monostatic  configuration, 
returns  are  localized,  occurring  at  the  same  time  position  (range)  from  frame-to-frame. 
Removal  is  relatively  easy  with  the  use  of  ensemble  averaging. 
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Here  we  consider  the  collection  of  a  fixed  number  of  data  frames,  and  perform  a  running 
average  to  reduce  noise  fluctuations.  For  each  block  of  M  frames,  the  resultant  is  then 
subtracted  from  successive  data  sets  to  eliminate  stationary  signal  components.  Note  that  for  a 
bistatic  or  multistatic  configuration,  application  of  this  method  is  not  as  simple  because  the 
relative  positions  of  scatterers  will  change  depending  on  the  particular  Tx  or  Rx  antenna 
element. 

As  an  example  of  how  fluctuations  are  reduced,  we  consider  a  white  noise  source,  as 
portrayed  in  Figure  8.  The  effects  of  the  averaging  are  summarized  in  Table  1.  In  this  context, 
N  denotes  the  number  of  frames  averaged  to  produce  the  trace.  A  running  average  is  obtained 
with  the  following  lines  of  Matlab  code: 

for  k  =  1  :N  %  loop  over  #  averages 

X  =  rand(2000, 1 ):  %  random  vector 

sum  =  sum  +  X; 

av  =  sum/k;  %  compute  running  average 

end 


This  “ideal”  white  noise  source  acts  as  a  reference  for  comparison  with  collected  radar  data. 


N  =  1 


N  =  10 


N  =  100 
*  ********»' 


1 


N  =  1000 


Figure  8.  Example  of  averaging  to  improve  SNR.  The  plots  show  how  averaging  reduces  noise  variance 

to  extract  a  signal  of  interest  (units  are  arbitrary). 

In  this  example,  we  see  that  an  arbitrary  signal  of  constant  amplitude  is  initially  undetectable 
given  the  noise  background.  As  the  number  of  frames  averaged,  A,  increases,  the  signal  to 
noise  ratio  (SNR)  continually  improves.  By  reducing  noise  variance,  the  averaging  provides 
an  effective  gain  of  approximately  G  «  101og10  N ,  with  respect  to  a  single  pulse. 
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Table  1.  Representative  processing  gains  achieved  from  noise  averaging. 


#  averages,  N 

noise  variance,  a2  [dB] 

SNR  [dB] 

1 

-4.8 

-4.3 

10 

-14.8 

5.7 

20 

-18.1 

9.1 

50 

-21.8 

12.7 

100 

-24.8 

15.8 

1000 

-34.6 

25.5 

3.5  SNR  improvement 

The  improvement  in  signal  quality  can  be  quantified  by  examining  the  signal-to-noise  ratio 
(SNR)  and  noise  variance.  The  improvement  factor  for  an  idealized  noise  source  has  already 
been  examined.  In  this  section,  we  consider  two  cases,  one  in  which  the  radar  data  is 
processed  by  bandpass  and  notch  filtering  and  is  averaged  over  a  number  of  frames  to  reduce 
fluctuations.  The  second  case  is  where  the  data  is  once  again  averaged,  but  is  not  filtered. 

The  variance  cr2 ,  follows  a  1/N  trend  as  is  expected  for  noise  power  reduction,  but  starts  to 
flatten  out  at  higher  values  of  N  (Figure  9). 


Figure  9.  Idealized  and  actual  variance  (signal  performance). 


This  deviation  can  be  explained  by  understanding  that  the  actual  returned  signal  is  not 
composed  entirely  of  thermal  fluctuations-there  are  components  from  the  target  that  have 
relatively  large  amplitude  in  comparison  to  the  background.  As  the  number  of  frames 
averaged  increases,  the  influence  of  the  target  returns  increases  with  respect  to  the  “pure” 


12 


DRDC  Ottawa  TR  2004-251 


noise,  resulting  in  a  greater  variance.  This  is  the  trend  displayed  in  Figure  9.  The  two  plots 
have  been  scaled  to  align  for  N=1 ,  a  single  pulse. 

An  idealized  SNR  can  be  calculated  as  SNR  =  -101og10(cr2),  where  the  reference  signal  is 

set  to  unity.  However,  actual  SNR  performance  (Figure  10),  shows  a  drastic  departure  from 
the  ideal  noise  source. 


Figure  10.  Idealized  and  actual  SNR  performance. 


For  the  radar  data  displayed  in  this  figure,  averaging  over  a  relatively  small  number  of  pulses 
is  insufficient  to  raise  the  desired  returns  appreciably  above  the  noise  floor.  Typical  SNR 
values  at  N  =  1  are  around  18  to  19  dB ,  and  the  trend  does  not  initially  follow  that  of  the  ideal 
noise  source.  A  possible  explanation  for  this  is  that  not  enough  signal  energy  is  collected  for 
lower  values  of  N.  The  idealized  curve  is  white  noise,  something  only  approached  for  large  N 
in  the  case  of  collected  data;  the  overall  result  is  a  smaller  target  signal  response  in  the 
presence  of  relatively  large  noise  contributions. 

Beyond  this,  SNR  improves,  and  does  so  at  a  rate  comparable  to  white  noise.  In  this  particular 
instance,  SNR  peaks  at  approximately  35  dB ,  an  overall  improvement  of  -15  dB.  For  greater 
values  (for  example,  N  >  80),  the  averaging  takes  place  over  a  long  enough  time  period  that 
energy  is  spread  due  to  target  motion.  This  effectively  defines  the  upper  limit  of  N  that  can  be 
achieved  without  additional  processing  to  focus  the  returns.  In  Figure  10  this  is  the  point 
where  the  calculated  SNR  flattens  out,  and  then  begins  to  decrease  with  N.  The  effects  on  the 
pulse  itself  are  observed  in  Figure  1 1 . 
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N  =  10 


unfiltered 
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filtered 


N  =  500 


v^M(V^v 


Figure  11.  Received  signal  traces  (N=10,  100,  500;  range  =  21m). 

This  diagram  contains  information  on  a  variety  of  trends.  There  are  three  pairs  of  plots,  with 
pulse  returns  for  TV  =  10,  100  and  500.  In  turn,  each  pair  shows  the  outcome  of  unfiltered 
processing  (top)  and  filtered  processing  (bottom).  That  is  .filtering  includes  the  application  of 
a  bandpass  filter  for  1  to  3GHz,  and  background  noise  conditioning  as  described  previously. 
Examining  the  trend,  it  is  obvious  that  filtering  eliminates  much  of  the  underlying 
interference.  For  example,  a  low-frequency  oscillation  is  plainly  visible  in  the  topmost  trace, 
but  the  combination  of  bandpass  and  notch  filtering  removes  this.  To  remove  the  higher 
frequency  fluctuations,  we  employ  the  ensemble  averaging  technique,  with  the  outcome 
displayed  in  the  bottom  traces.  As  N  increases,  more  and  more  of  the  fluctuations  are 
removed,  but  the  target  return  begins  to  disperse  again,  due  to  target  motion.  The  trace  for  N  = 
500  corresponds  to  the  rightmost  point  of  Figure  10,  where  the  calculated  SNR  is  far  past  its 
peak  value. 

A  more  quantitative  description  of  this  is  contained  in  Figure  12,  comparing  SNR  and 
variance  as  functions  of  N  for  both  filtered  and  unfiltered  data.  From  the  processing  described 
here,  we  are  able  to  improve  the  signal  to  noise  ratio  from  a  minimum  of  ~  3dB  at  N  =  1  to  a 
maximum  of  7.4dB  at  N=  40. 
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variance,  a2  [dB] 


Figure  12.  Variance  and  SNR  vs.  N  for  both  filtered  and  unfiltered  data. 


Figure  13  contains  an  example  of  improved  resolution,  showing  the  fine  structure  of  a  pulse 
return.  As  much  of  the  processing  involves  analytic  signals,  the  amplitude  (complex  envelope) 
is  superimposed  on  the  trace. 


Figure  13.  Comparison  of  return  pulse  waveform  and  its  complex  envelope  (via  Hilbert  transform). 
Given  a  real-valued  signal  waveform  s(t ),  the  corresponding  analytic  signal  is  defined  as 

i//(t)=s(t)  +  iH[s(t)],  (2) 

where  the  imaginary  component,  H[s(t )]  ,  is  the  Hilbert  transform  of  s(t) — 
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H[s(t)\  =  <P 


(3) 


I  r  ^ %t > 

n 

As  this  is  an  improper  integral  by  virtue  of  the  singularity  at  t  =  t',  it  is  expressed  in  terms  of 
a  Cauchy  principle  value,  denoted  by  the  symbol  <P : 


(p 


1 r 42*1 

=  —  lim 

r^f'+r  s(,Vl 

n  £->° 
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>Ki 

+ 

Vi 
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>l-*i 
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I 

(4) 


This  transform  remains  entirely  within  the  time  domain.  Further  details  regarding  the  use  of 
Hilbert  transforms  in  SAR  image  processing  are  provided  in  [1]. 

Figure  14  compares  returns  from  the  same  target  (rotating  sphere)  at  different  distances 
ranging  from  14  to  49  meters. 


0  1  2  3  4m 

Figure  14.  Complex  envelope  (Hilbert  transform)  for  pulse  returns  at  various  target  ranges. 
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4.  The  impulse  radar  signal 


In  contrast  with  narrowband  signals,  UWB-SP  waveforms  are  baseband.  They  are  created 
using  impulse  generators,  and  have  no  explicit  carrier  frequency.  To  characterize  returns,  we 
can  refer  to  center  frequencies  and  bandwidths,  but  keep  in  mind  these  are  controlled  only 
indirectly  by  variation  of  pulse  width  and  rise/fall  times,  which  are  in  turn,  hardware 
dependent. 

4.1  Gaussian  model 

To  simplify  the  analysis,  we  consider  the  initial  impulse  to  be  in  the  form  of  a  Gaussian  pulse. 
This  is  a  reasonable  assumption  based  on  comparison  with  a  typical  pulse  obtained  directly 
from  the  generator4  (Figure  15) 


Figure  15.  Measured  data  (solid  line)  overlaid  with  Gaussian  pulse  (dashed  line). 
FWHM  pulse  width  is  r  =  0.42ns. 


4  Avtech  Model  AVH-HV1-C-T1.  This  model  can  output  a  variable  level  pulse  up  to  100  V  peak  with  a 
width  of  600ps.  Voltage  was  attenuated  to  keep  levels  within  the  limits  of  the  oscilloscope.  The  actual 
output,  taking  into  account  total  attenuation  and  amplitude  setting,  was  8.6  V. 
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For  this,  we  have  performed  a  fit  on  the  data,  based  on  a  Gaussian  distribution 


1  («-* o)2 

At)  =  —fT=e  2ff2  •  (5) 

An  easily  accessible  experimental  parameter  is  the  full-width  at  half-maximum  (FWHM) 
pulse  width.  In  the  following,  we  rewrite  the  Gaussian  pulse  in  terms  of  the  FWHM,  x. 

The  standard  deviation  crand  pulse  width  rare  related  by  the  relation 

r  =  2<j^2Ln2.  (  6  ) 

Introducing  the  parameters  A  =  (^t42jt  )  and  a  =  (2 y2r2 )  ,  along  with  the  scaling 
y  =  *  0.42  ,  we  obtain 


f(t)=Ae~a{t~tof .  (7) 

We  find  that  for  the  measured  data  in  Figure  15,  the  FWHM  pulse  width  is  0.42ns.  This 
appears  to  be  an  extremely  narrow  pulse,  but  recall  this  was  collected  directly  from  the  pulse 
generator,  and  experienced  none  of  the  distorting  factors  typical  of  the  transmit  path. 

In  choosing  a  Gaussian  form,  we  lose  some  flexibility  for  system  parameters,  as  we  cannot 
independently  specify  pulse  width  and  rise  times.  In  narrowband  pulsed  radar,  rise  time  is  not 
typically  considered  and  signal  bandwidth  is  defined  as  1/At ,  where  At  is  the  pulse  width.  As 
the  pulse  width  decreases  however,  the  rise  time  becomes  a  significant  fraction  of  the  total 
pulse  width.  As  rise  time  decreases,  the  pulse  must  comprise  higher  frequency  components. 

We  can  however,  use  this  model  to  examine  pulse  returns  and  to  estimate  transmitted  pulse 
widths.  A  previous  paper  [1]  established  that  target  return  waveforms  could  be  modeled  as 
Gaussian  fourth-order  derivatives. 

Antenna  calibration  data  is  supplied  by  the  manufacturer,  Farr  Research  Inc.  [9],  in  the  form 
of  an  impulse  response  h(t ),  where  Vrec(t)  =  h(t)  *  Eincident(t).  Starting  with  the  Gaussian 
waveform  as  the  incident  pulse,  this  is  convolved  with  the  impulse  response  of  the  transmit 
antenna, 


PtA)  =  P*K  =\p{r)hTx{t-r)dr ,  (1) 

with  the  result  shown  in  Figure  16a.  In  turn,  this  result  undergoes  a  180°  shift  to  model  a 
simple  reflection  from  a  scatterer,  and  is  finally  convolved  with  the  receive  antenna  impulse 
response  to  yield  Figure  16b. 


pM=  ~Ptx  *hRx=~\ Ptx  W hRx{t-r)dT, 


(2) 
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Figure  16.  Predicted  Tx  and  Rx  signals  resulting  from  convolution  analysis. 

This  result  can  be  compared  with  the  target  response  from  measured  data  as  displayed  in 
Figure  17.  There  is  obviously  qualitative  agreement,  although  we  would  not  expect  an  exact 
match  owing  to  characteristics  of  the  physical  environment  and  the  specific  target  RCS. 

0.03 

0.02 

0.01 

2 

3  o 

-0.01 

-0.02 


distance  [m] 


Figure  17.  Target  response  (measured  data) 


Referring  back  to  the  expression  for  f(t ),  and  obtaining  the  fourth-order  derivative, 

f(4](t)  =  A(l6a4t4  -48a3t2  +  I2a2)e~°“2 , 
from  which  we  can  derive,  via  the  Fourier  transform, 

F(4Hv)  =  167r4v4A^e-’r2y2/a, 

V  a 


where  the  frequency  in  Hz  is  denoted  v. 


(8) 


(9) 
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4.2  Model  validation 


In  the  following,  we  compare  this  model  with  experiment.  Figure  18  shows  the  time-domain 
pulse  and  how  it  reproduces  the  shape  and  duration  of  experimental  data.  The  only  free 
parameter  is  r,  the  FWHM  pulse  width  of  the  original  Gaussian  pulse  as  in  Figure  15. 


Figure  18.  Gaussian  4th  order  time-domain  pulse  (dotted  line)  overlaid  on  experimental 

data  (t=  1.3  ns  FWHM). 

Based  on  this  model,  we  can  obtain  an  analytic  value  for  the  power  spectrum.  Once  again,  the 
agreement  with  experiment  is  quite  good.  The  pulse  broadening  that  occurs  between  the 
original  pulse  generator  output,  x  =  0.42  ns  (Figure  15)  and  the  received  pulse  (FWHM  fitted 
to  1.3  ns,  Figure  18)  is  presumably  due  to  bandwidth  limitations  of  the  RF  components  and 
sampling.  The  same  effect,  with  comparable  magnitude  has  been  noted  elsewhere  ( [10]  and 
[11]): 

•  x  =  0.66  ns  broadened  to  2.5  ns; 

•  x  =  0.66  ns  broadened  to  1 .44  ns. 
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As  a  further  comparison  between  theory  and  experiment,  Figure  19  shows  the  power  spectra 
corresponding  to  the  time-domain  pulse  in  Figure  18.  The  red  dotted  line  again  refers  to 
theoretical  prediction.  The  match  for  frequencies  under  1GHz  is  very  good,  providing  more 
support  for  the  waveform  interpretation  elaborated  on  in  this  report.  The  spectrum  beyond 
1GHz  is  predominantly  noise  generated. 


Center  frequency  fc  =  0.38GHz,  FWHM3dB  =  0.22GHz. 
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4.3  Radiation  patterns 


The  standard  approach  to  analyzing  radiation  patterns  is  based  on  transmission  of  sinusoidal 
waveforms.  This  is  no  longer  the  case  for  UWB  impulse  radar:  with  no  carrier  and  a  narrow 
pulse  width,  the  signal  is  distorted  as  a  consequence  of  transmission.  This  becomes  an 
important  consideration  in  analysis  of  waveforms  [12]. 

The  electric  field  intensity  at  any  point  is  E  =  -VV  -  8 tA  ,  where  V  is  the  electric  scalar 

potential  and  A  the  magnetic  vector  potential,  defined  by  B  =  V  x  A  (see  for  example,  [13]). 
Next,  take  the  time  derivative 


8B_ 

8t 


f  8A^ 


\dtj 


This  is  substituted  into  Faraday’s  law  (MKS) 


VxE  =  - 


dB_ 

dt 


-Vx 


dt 


(10) 


(11) 


which  implies  E  =  8 tA  .  In  turn,  the  magnetic  vector  potential  is  proportional  to  current. 

From  these  results,  we  see  that  the  field  components  from  an  electric  dipole  are  proportional 
to  the  time-derivative  of  the  current  l{t  \  and  to  the  second  time-derivative  of  the  input 

voltage  v{t  \  The  notation  t  —t  —  r  /  c  indicates  retarded  time.  Since  the  operators  are 
linear,  the  combined  operations  of  transmit  and  receive  will  yield  the  fourth  time-derivative  as 
obtained  with  the  Gaussian  model. 
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5.  System characterization 


The  basis  of  the  present  study  is  the  radar  equation,  which  relates  the  received  power  to 
various  parameters  associated  with  EM  propagation  and  detection.  The  following  analysis 
uses  a  realistic  parameter  set,  although  the  calculations  should  be  interpreted  carefully.  In 
specific  scenarios,  it  may  be  possible  to  enhance  the  radar  performance  beyond  the 
contemplated  model  [5]. 

The  radar  equation  below  provides  the  SNR  for  multiple  pulses  of  power,  Ph  as  a  function  of 
antenna  gain,  G,  wavelength,  X  (nominally  at  the  center  frequency),  target  RCS,  cr0 ,  range,  7?, 

system  losses,  Z,  noise  figure,  Fm  and  bandwidth,  B.  Boltzmann’s  constant5  is  denoted  by  k 
and  T0  is  the  reference  temperature,  290K.  Degradation  due  to  clutter  returns  is  captured  with 
the  parameter  Lc.  The  value  N  denotes  the  number  of  pulses  integrated  to  obtain  the  given 
SNR. 


With  these  parameters,  the  signal  to  noise  ratio  becomes  [14] 


SNR  = 


PtG2A2cr0N 

(4^y¥lfJPblc  ‘ 


(12) 


The  theoretical  antenna  gain  is: 


G  = 


A2 


(13) 


where  Ae  =  rjA  is  the  effective  antenna  area.  The  antenna  efficiency  is  characterized  by  the 
parameter  r|,  such  that  0  <  77  <  1 . 

The  radar  power  is  derived  from  a  pulse  voltage  of  ~103V  (peak)  at  the  terminals  of  a  line  of 
characteristic  impedance  of  between  50  and  100  Q.  The  antenna  gains  are  based  on  a  dish 
diameter  of  about  0.5.  The  estimated  losses  comprise  10  dB  in  the  antenna  feed  and  the 
wideband  antenna  itself  and  occur  on  both  transmit  and  receive.  The  high  bandwidth  is 
obtained  by  heavy  damping  designed  to  broaden  antenna  response  without  generating  ringing 
or  phase  dispersion. 

The  noise  figure  includes  the  effect  of  noise  introduced  by  the  front-end  amplifier.  The  RCS  is 
that  of  a  0.15  m  radius  sphere  {ttt2  «  0.07m2).  In  practice  the  returns  may  be  distributed 
among  many  resolution  cells;  for  example,  if  the  return  is  spread  over  10  cells,  the  effective 
RCS,  aeff,  for  any  one  cell  will  be  <Jeff  «  0.  lcr0 .  The  pulse-integrated  SNR  predicted  by  the 

radar  equation  is  given  in  the  last  row  of  the  table,  along  with  the  value  derived  directly  from 
measured  data.  Table  2  provides  a  summary  of  the  parameter  set. 


5  k  =  1.38xl0"23  J/K,  or  -228.6  dBw/Hz/K. 
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Table  2.  Radar  equation  parameters 


Parameter 

Value 

Symbol 

Notes  and  information  source 

Peak  transmit  power  [w] 

104 

Pt 

Instrument  spec 

Center  frequency  [GHz] 

0.4 

fc 

Experiment  /  data  fitting  (see  page  21 ) 

Bandwidth  [GHz] 

2 

B 

Experiment  /  data  fitting  (see  page  21) 

Wavelength  [m] 

0.75 

A, 

Calculated  from  fc 

Antenna  gain  [dBi] 

14.8 

G 

Experiment  (see  page  27) 

Antenna  efficiency 

0.5 

n 

Estimated 

(typical  -  antenna  resistively  loaded) 

Antenna  area  [m2] 

0.19 

A 

Measured 

RCS  [m2] 

0.07 

ft  r2 sphere  (^sphere  ~  0.15m) 

Target  is  metallic  sphere 

Range  [m] 

28 

R 

Measured 

System  losses  [dB] 

16 

L 

Estimated  (cable,  connector  losses) 

Clutter  [dB] 

3 

Lc 

Estimated  (monostatic  configuration) 

Receiver  noise  figure  [dB] 

10 

Fn 

Instrument  spec  /  estimate  (see  page  28) 

Pulse  integration 

50 

N 

Processing  parameter 

SNR  (calculated)  [dB] 

43.7 

SNR  (measurement)  [dB] 

44.0 

This  parameter  set  was  used  to  construct  a  set  of  SNR  vs.  Range  curves  which  are  in  turn, 
plotted  along  with  a  set  of  discrete  data  points  derived  from  experiment  (refer  to  Figure  20). 
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Figure  20.  Predicted  SNR  vs.  Range  (solid  &  dashed  curves)  compared 
with  experimental  results  (points). 

There  are  two  curves  displayed,  for  different  antenna  gains.  Because  the  antenna  gain  appears 
as  a  squared  quantity  in  the  range  equation,  it  is  a  relatively  sensitive  parameter  as  can  be 
observed  in  the  plots  -  the  predicted  SNR  varies  across  a  range  of  4dB  for  AG  =  2dB.  The 
method  for  estimating  antenna  gain  is  outlined  in  the  next  section. 

The  apparently  poor  performance  of  the  model  at  a  shorter  range  can  be  misleading — it  is 
actually  the  result  of  the  transmit  pulse  not  being  adequately  time-gated  out,  and  excessive 
power  gets  into  the  receive  antenna.  Although  the  signal  strength  is  good,  there  is  excessive 
clutter,  raising  the  noise  floor.  This  in  turn,  artificially  lowers  the  signal  to  noise  ratio. 

5.1  Antenna  beamwidth  and  gain 

To  obtain  an  estimate  of  the  Farr  antenna  beamwidth,  a  set  of  measurements  were  taken  to 
record  peak  amplitude  over  a  range  of  aspect  angles.  In  what  follows,  aspect  angle  is  defined 
as  that  angle  existing  between  the  main  beam  of  the  transmitting  antenna  and  the  main  beam 
of  the  receive  antenna.  That  is,  for  one  antenna  pointed  directly  at  the  other,  the  aspect  angle 
is  zero  degrees. 
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Transmit-Receive  (TX/Rx)  separation  remained  fixed  for  each  measurement,  and  the  Tx 
antenna  was  manually  swept  through  a  range  of  aspect  angles  from  0°  to  60°. 


Table  3.  Data  for  received  signal  strength  vs.  aspect  angle. 


Aspect  angle 
[degrees] 

Receive  voltage  [mV] 

(on  channel  2) 

Gain  [dB] 

0 

522.6 

0.0 

10 

423.7 

-1.82 

20 

291.3 

-5.08 

30 

242.2 

-6.68 

40 

206.4 

-8.07 

50 

181.6 

-9.18 

60 

159.2 

- 10.32 

Transmit  voltage  recorded  on  channel  1  was  1.39  V.  Assuming  a  unit  resistance,  we  can 
calculate  the  received  power  at  each  angle,  P{0)  =  m  V x  IQ-3/!!  .  From  this,  the  relative  gain 
is 


Gj6>)=  101oglo 


f 


V 


max  J 


(14) 


In  Figure  21,  these  values  were  fitted  to  a  quadratic  function,  yielding  an  approximate  value 
of  26°  for  the  3dB  point  of  the  antenna  beamwidth. 


[degrees] 


Figure  21.  Antenna  gain  vs.  aspect  angle. 
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Another  way  of  viewing  this  information  is  on  a  polar  plot,  (refer  to  Figure  22)  where  the 
beam  pattern  is  more  apparent. 


90 


Figure  22.  Polar  plot  showing  experimental  antenna  beamwidth. 

We  can  use  the  measured  beamwidth  to  obtain  an  estimate  of  antenna  gain.  To  do  so,  we  first 
consider  the  surface  area  of  a  sphere, 


In  n 

sin  6  dO  =  An  .  (15) 

surface  0  0 

Converting  from  radian,  this  is  4 n  x  (l  80  /  7r)2  =  41,253  square  degrees.  An  antenna  with  3- 
dB  horizontal  and  vertical  beamwidths  of  0H  and  6V  respectively,  therefore  radiates  into 
6 H0V  square  degrees  out  of  the  total  surface  area.  Accounting  for  an  antenna  efficiency  of  rj , 
we  have 


G~ 


41,253 

- 77 

Wv 


(16) 


For  example,  given  the  symmetry  of  the  Farr  antenna  under  consideration,  we  expect  the 
horizontal  and  vertical  beamwidths  to  be  equivalent.  For  the  moment,  we  neglect  ground 

effects.  Therefore,  0H  =  6V  —  26° ,  and  setting  r/  =  0.5  results  in  the  approximate  value 

G  «  (41,253/  262)o.5  =30.5,  or  in  dB,  GdB  «  14.8 dBi.  Here,  dBi  refers  to  gain  with  respect 
to  isotropic ,  an  antenna  that  radiates  equally  in  all  directions. 
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Due  to  the  experimental  method,  there  exists  some  uncertainty  in  the  beamwidth;  the 
estimated  gain  should  be  considered  an  upper  limit  until  more  precise  measurements  can  be 
made. 


5.2  Receiver  noise  figure 

In  the  series  of  measurements  carried  out  here,  the  Tektronix  oscilloscope  effectively  acts  as  a 
receiver  and  detector.  As  such,  it  is  subject  to  noise  impairments,  which  we  discuss  presently. 

From  the  published  specifications  of  the  TDS  7404  [15],  the  highest  sensitivity  is 
2mV/division,  resulting  in  16mV  at  full-scale  deflection.  In  turn,  the  hardware  employs  an 
eight-bit  digitizer.  At  a  sampling  rate  of  20GS/s  and  analog  bandwidth  of  4GHz,  the  A/D 
resolution  of  the  scope  is  limited  to  4.8  effective  bits.  This  means  that  approximately  3  bits  are 
lost  due  to  quantization  noise,  jitter,  A/D  nonlinearity,  and  interpolation  error. 

For  the  full  8-bits  (28  =  256  levels),  we  have 

16  mV  /  256  levels  =  0.0625  mV/level. 

Allowing  3  bits  for  noise,  (23  =  8  levels), 

8  *  0.0625  mV/level  =  0.5mV. 

In  the  calculation  of  noise  figure  this  becomes  the  output  voltage.  The  input  signal  is  given  by 
thermal  noise  [16]: 

V„l=4  kTBR 

=  4  (l  .38  x  1 023  J  /  /O(300ii:)(4.0  x  1 09  //z)(50Q) 

=  3.3x10'9F2 


Therefore,  the  oscilloscope  noise  figure  (NF)  is 


Q.25otF2 
3.3x1 0~9mV2 


75.8  [NFdB  « \9dB). 


5.3  SNR  and  path  loss 

Figure  23  shows  some  aspects  of  the  transient  behaviour.  The  upper  plot  is  a  trace  of  target 
range  vs.  scan  (pulse)  number  -  that  is,  for  each  pulse  transmitted  and  received,  the  echoed 
waveform  is  analyzed  to  determine  the  range  associated  with  the  highest  return  voltage  peak. 
The  lower  trace  shows  how  the  maximum  amplitude  varies  with  scan  number. 

What  is  interesting  about  this  is  how  the  amplitude  drops  close  to  zero  for  the  points 
corresponding  to  extremes  of  the  range  sinusoid.  With  a  little  thought,  it  becomes  obvious  that 
this  is  due  to  the  fact  that  radial  velocity  vanishes  at  the  two  extrema  and  therefore  the  pulse- 
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to-pulse  subtraction  will  not  produce  an  appreciable  voltage  signal  (This  will  be  discussed 
further  in  Section  6,  Velocity  processing). 
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Figure  23.  Target  range  vs.  scan  number  (top),  Receive  signal  peak  vs.  scan  number  (bottom). 

However,  the  fluctuating  voltage  can  be  used  to  extract  signal  power,  and  in  turn,  estimate 
SNR.  The  lower  trace  of  Figure  23  was  averaged  to  obtain  a  mean  voltage  (Indicated  by  the 
horizontal  line  extending  across  the  trace).  A  more  accurate  value  could  be  obtained  by 
weighting  the  average  to  reduce  the  influence  of  the  two  voltage  drops,  but  this  was  not 
carried  out. 
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The  results  are  shown  in  Table  4,  where  the  SNR  over  a  set  of  target  ranges  was  determined. 
The  noise  power  was  measured  separately. 


Table  4.  Calculated  return  signal  power  and  SNR  (N  =  50;  noise  floor  =  -174  dBw/Hz) 


Range  [m] 

<Vrx>  [VOlts] 

Power  @  50Q  [dBw] 

SNR  [dB] 

14 

0.0248* 

-39.1 

49.9 

21 

0.0123* 

-45.2 

43.8 

28 

0.0224 

-50.0 

44.0 

35 

0.0155 

-53.2 

40.8 

42 

0.0097 

-57.3 

36.7 

49 

0.0078 

-59.1 

34.8 

*  A  10  dB  PAD  was  placed  in  the  receive  path  to  shield  the  oscilloscope.  Otherwise,  signal  power  received  at  these 
distances  could  damage  instrumentation.  To  get  proper  voltage,  scale  by  Vi  0  ~  3.2. 


As  a  means  of  verifying  these  results,  the  received  power  was  plotted  vs.  range  and  the  path 
loss  exponent  determined.  As  the  experimental  configuration  was  essentially  free-space 
propagation,  the  exponent  should  be  close  to  4,  as  shown  in  Figure  24.  The  dashed  line  shows 
strict  R~4  behaviour  as  comparison  (Normalized  for  14m  reference  range). 
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6.  Velocity  processing 


Traditionally,  the  processing  of  moving  targets  is  carried  out  in  the  frequency  domain  with 
data  analyzed  in  terms  of  the  Doppler  shift  in  frequencies: 


where  fD  is  the  Doppler  frequency  with  target  velocity  u.  The  wavelength  of  the  (narrowband) 
signal  is  denoted  X.  The  principle  of  operation  here  is  that  moving  targets  can  be  detected  in 
clutter  due  to  frequency  discrimination.  Stationary  targets  will  not  cause  a  frequency  shift  in 
the  carrier. 

6.1  Radar  imaging 

The  range  processor  described  in  previous  sections  was  modified,  yielding  a  2D  SAR 
processor.  Measurements  were  carried  out,  examining  the  motion  of  a  target  rotating  about  a 
radius  of  0.5m.  The  use  of  a  sphere  as  an  idealized  target  ensured  a  relatively  uniform  RCS, 
independent  of  the  point  of  its  orbit. 

Radar  images  for  targets  at  ranges  of  7,  21,  35  and  49m  are  displayed  in  Figure  25  and  Figure 
26.  Each  consists  of  the  reconstructed  image  and  its  associated  contour  diagram.  Using  the 
space-time  backprojection  technique  ([1]  and  [17]),  the  images  clearly  show  the  oscillatory 
motion  of  the  target,  but  the  signal  vanishes  as  radial  velocity  approaches  zero  at  the  two 
extremes  of  motion.  Note  that  motion  of  the  target  was  slow  enough  that  only  one  complete 
rotation  was  captured  for  imaging. 
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#  scans 


Figure  25.  ISAR  images  of  rotating  target  at  ranges  a).  7m  and  b).  21m. 
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Figure  26.  ISAR  images  of  rotating  target  at  ranges  a).  35m  and  b).  49m. 

A  comparison  of  contour  plots  for  a  target  at  21m  is  shown  in  Figure  27.  The  received  pulses 
have  been  averaged  from  20  to  200  data  frames.  The  most  obvious  trend  is  that  the  gap 
corresponding  to  near  and  far  extremes  tends  to  widen  as  the  averaging  is  reduced.  This  will 
be  discussed  further  in  section  6.2,  Non-Doppler  MTI. 
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N  =  20  N  =  50  N  =  100  N  =  200 


Figure  27.  Contour  plots  for  the  same  target  at  range  =  21m,  but  with  pulses  averaged  over  a  number  of 

values  of  N:  20,  50,  100,  and  200. 

6.2  Non-Doppler  MTI 

As  indicated  previously,  UWB-SP  radar  transmissions  are  carrier-free  and  we  are  unable  to 
exploit  Doppler  processing.  However,  target  motion  is  still  observable  by  examining  changes 
in  target  range  from  pulse  to  pulse.  The  basis  for  processing  now  becomes  the  change  in 
arrival  time  of  the  pulse.  In  effect,  we  employ  a  form  of  Area  Moving  Target  Indicator  (Area 
MTI),  comparing  the  complex  envelopes  of  successive  scans  to  identify  target  motion. 

We  consider  a  target  moving  at  velocity  u.  For  a  radar  transmitting  pulses  at  a  PRF  equal  to 
1/T,  the  target  will  move  a  distance  of  uT  for  each  transmitted  pulse.  The  ensuing  time  delay 
is  therefore  twice  the  distance  moved  by  the  target,  divided  by  the  propagation  velocity, 

8  =  2d/c  =  2 uT /c.  Furthermore,  when  each  pulse  is  transmitted,  an  echo  is  subsequently 
received  and  sampled  at  time  increments  of  At .  In  the  data  studied  here,  we  have 
At  =  5.0  x  10”11  v  (or  50ps).  Note  that  At  corresponds  to  a  range  bin. 

To  obtain  the  detectable  velocity,  we  require  the  target  to  move  at  least  k  range  bins  during  the 
observation  time  of  N  pulses  [18]: 
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•  For  velocity  u ,  total  distance  is  NTu  . 

•  For  k  range  bins,  the  total  distance  is  k{cAtj 2) 
Equating  the  two,  we  have 


u  = 


kcAt 

2NT 


(18) 


To  obtain  the  minimum  detectable  velocity,  we  require  the  target  to  traverse  at  least  two  range 
bins  (k=2).  Otherwise,  the  target  is  considered  stationary. 


cAt  _  cAt 
~NT  ~  "a” 


PRF 


(19) 


The  experimental  setup  consists  of  the  target  centered  at  R0,  rotating  with  angular  velocity  c o. 
The  radius  of  the  motion  is  r.  As  the  target  traces  out  a  circular  path,  the  maximum  range  from 
the  radar  is  R0  +  r  at  cp  —  0,  while  the  minimum  range  is  R0—r  at  (p  =  n  (note  that  (p  =  cot). 

In  general,  the  range  is  described  by  the  position  vector  R  =  R0  +  r  ,  as  indicated  in  Figure  28. 


Figure  28.  Configuration  for  measurement  of  oscillatory  target  motion. 

(co  =  1.2  rpm  ) 

Using  the  law  of  cosines,  we  can  express  the  radial  distance  from  antenna  to  target,  7?,  as  a 
function  of  cp. 
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(20) 


R  =  /r^  +  V2  -  2 i?0r  cos(^  +  ;r) 

=  ^/i?02  +  r2  +  27?0  r  cos(y>) 


The  radial  velocity,  in  turn,  is  given  by 

dR 


2coR0rsm{(p) 


1  radial 


dt  2-JrI  +r2  +2i?0rcos(^) 
The  maximum  velocity  occurs  at  the  points  cp  =  ±;r/2,  and  we  have 


max 
4  radial 


dR 

a>R0r 

dt 

A2 +}~2 

(21) 


(22) 


This  expression  is  further  simplified  by  making  the  reasonable  assumption  that  i?0  »  r, 
which  results  in  ~  cor- 

To  make  contact  with  experimental  data,  the  target  executes  one  rotation  about  a  0.5m  radius 
every  50s.  This  is  equivalent  to  an  angular  velocity  of  1.2  rpm,  or  co  =  0.125  rad/s.  The 
maximum  radial  velocity  is  therefore  0.0 63  m/s.  Once  we  have  scaled  the  displacement  (as 
indicated  in  Figure  29),  we  can  determine  the  corresponding  velocity  (Figure  30). 


/=  20.97  +  0.5*sin(0.065*x  -19.7) 


Figure  29.  Experimental  data  and  parametric  fit  used  to  extract  minimum  velocity  threshold. 
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Figure  30.  Displacement  end  points  mapped  onto  radial  velocity  function. 

In  this  particular  example,  the  PRF  is  100Hz,  and  the  number  of  pulses  integrated,  N  =  50. 
From  the  expressions  developed  previously,  we  have  for  the  minimum  detectable  velocity  — 

vmm=cAtPRF/N 

(2.998 xl08m/4  (5.0xl0“ns)  (lOOHz)  n  Mn  , 

=  ± - ^ 0.029m/ s 

50 

Figure  31  compares  theoretical  values  with  those  obtained  from  processed  images.  While  the 
trend  follows  theory,  there  is  a  discrepancy  at  lower  pulse  averaging.  This  is  presumably  a 
result  of  the  increased  variance  for  low  averaging  and  the  resulting  error  in  determining  the 
threshold  value. 


Figure  31.  Minimum  detectable  target  velocity  as  a  function  of  PRF  and  #  averages. 
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7.  Wall  penetration  losses 


As  concrete  was  selected  to  be  the  building  material  of  interest  for  this  project,  we  need  to 
determine  the  effects  of  moisture  content  on  signal  propagation  and  attenuation. 

7.1  Wall  construction 

The  materials  used  in  this  experiment  were  standard  concrete  blocks  of  approximate 
dimensions  20  x  20  x  40  cm,  with  two  air  pockets  inside.  Recent  studies,  [19]  and  [20],  report 
that  dry,  solid  concrete  exhibits  a  dielectric  constant  of  8  to  12.  The  presence  of  air  gaps  in  the 
blocks  will  tend  to  reduce  that  value.  However,  to  our  knowledge,  little  has  been  mentioned  in 
the  literature  regarding  moisture  content  and  to  what  extent  it  will  increase  attenuation. 

A  concrete  block  wall  as  depicted  in  Figure  32  was  constructed  to  obtain  transmission  loss 
data.  The  transmit  antenna  was  moved  from  2  ft.  to  5  ft.  (60  to  152cm)  in  two-inch  increments 
(5cm),  while  the  receive  antenna  remained  fixed  at  a  distance  of  2  ft.  (60cm)  from  the  wall. 
From  this,  we  were  able  to  obtain  a  set  of  closely  spaced  transmission  data  over  a  range  of 
approximately  six  wavelengths. 


T  =  19  cm 
H  W  H - 


Figure  32.  Receive/Transmit  configuration  for  measurements  (side  view). 


The  wall  dimensions  are  53”  x  74.75”  x  7.5”  (135  x  190  x  19  cm).  With  an  antenna  height  of 
98  cm  and  beamwidth  of  approximately  26°,  there  should  be  minimal  direct  signal  power 
around  the  wall  and  limited  diffraction  phenomena.  In  fact,  the  transmit  antenna  would  need 
to  be  moved  about  nine  feet  from  the  wall  before  leakage  of  this  sort  became  a  factor. 
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7.2  Data  collection 


Following  this,  three  sets  of  measurements  were  taken: 

1 .  dry /damp  bricks  —  Bricks  were  stored  outside  subject  to  both  ambient 
temperature  and  climate.  The  bricks  were  not  visibly  wet,  but  will  exhibit  typical 
moisture  content.  No  attempt  was  made  to  prepare  a  reference  set  of  oven-baked 
dry  bricks. 

2.  wet  bricks  —  Bricks  were  thoroughly  soaked  by  repeatedly  hosing  the  wall  down. 
Water  was  applied  before  each  measurement  to  ensure  uniform  results.  This 
should  provide  an  upper  bound  on  the  degree  of  attenuation  due  to  water. 

3.  line  of  sight  (LOS)  —  The  wall  was  removed,  and  LOS  “free-space” 
measurements  taken  to  provide  a  reference. 

A  diagram  of  the  experimental  configuration  appears  in  Figure  33. 


T  DS  740*  Oscilloscope 


Figure  33.  Experimental  configuration 

Dividing  through-wall  values  with  the  corresponding  LOS  values  produces  normalized 
measurements,  eliminating  the  effect  of  1/r2  spreading  losses,  leaving  only  wall  attenuation 
loss.  A  typical  set  of  raw  data  is  shown  in  Figure  34. 
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time  [ns] 

Figure  34.  Unprocessed  data  from  wall  attenuation  measurements  (@  1.361). 

Tabulating  the  data  as  collected  in  Figure  34,  recording  only  peak  values  as  a  function  of 
distance  from  the  wall,  the  path  loss  due  to  line  of  sight  (LOS),  and  for  both  dry  and  wet 
cindercrete  blocks  is  depicted  graphically  in  Figure  35 


Figure  35.  One-way  propagation  loss,  unnormalized. 

Figure  36  shows  values  corrected  for  free-space  path  loss.  These  results  show  that  over  a 
range  of  about  six  wavelengths  the  through-wall  losses  are  between  7  to  9  dB  for  ambient 
blocks  and  between  9  to  1 1  dB  for  wet  blocks. 
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Figure  36.  Corrected  losses  for  wall  attenuation  (one-way). 

There  is  a  slight  upward  trend  in  the  data  that  seemingly  indicates  wall  attenuation  decreases 
with  distance.  A  possible  explanation  for  this  is  that  as  the  transmit  antenna  is  moved  farther 
away  from  the  wall,  a  portion  of  the  energy  “leaks”  around,  either  by  diffraction  or  some 
multipath  configuration  that  bypasses  the  attenuation  of  the  wall  material  itself. 

Note  that  this  is  a  systematic  effect,  which  can  be  verified  by  considering  the  differential  loss 
between  ambient  and  wet  blocks  (see  Figure  37)  is  relatively  flat.  The  average  loss  calculated 
from  this  is  Aloss  =  -2.16  +  0. 19  dB ,  with  a  slope  of  0.07  ±  0.02  dB/X. 
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7.3  Wall  reflections  and  waveform  structure 


Part  of  the  loss  incurred  when  transmitting  a  signal  through  a  wall  is  due  to  reflections  at  the 
air-dielectric  interface.  The  ratio  of  transmitted  and  reflected  signal  energy  is  given  by  the 
Fresnel  equations  [21].  These  can  be  expressed  in  quite  general  terms  for  both  electric  and 
magnetic  fields,  and  for  perpendicular  and  parallel  field  orientation.  For  the  purposes  of  this 
report,  it  is  sufficient  to  consider  the  more  restricted  case  of  normal  incidence. 

We  consider  perfect  dielectrics,  where  conductivity  is  neglected  ( o  «  0 );  all  materials  are 
non-magnetic  so  that  the  permeability  is  unity  (//  =  1).  Under  these  conditions  the  refractive 

index  takes  the  form  nf  «  ,  where  the  subscript  i  denotes  either  an  air  or  dielectric 

medium. 


It  is  convenient  to  make  calculations  in  terms  of  power,  where 


P  = 


nx  +  n2 


(23) 


is  the  fraction  of  reflected  incident  power.  The  transmission  coefficient,  from  conservation  of 
energy,  is  r  =  1  -  p  .  Power  losses  are  often  expressed  in  dB: 


Lp  =201og10 


n\  ~n2 
nx  +  n2 


and 


4  =201ogi0 


2k2 
nx  +  n2 


(24) 


(25) 


Consider  Figure  38,  depicting  the  one-way  transmission  of  a  signal  through  a  wall  with 
relative  dielectric  constant  £2 .  The  relative  dielectric  constant  of  air,  sx ,  is  approximately  one. 


—  For  an  incident  signal  of  power  Pt ,  the  power  transmitted  across  the  first  interface  into  the 
medium  is  Pm  =  rPt  (watts). 

—  In  terms  of  dB ,  this  is  Pm  dBm  =  PidBm  +  Lt  dB ,  where  dBm  is  a  convenient  measure  of  power 
with  respect  to  one  milliwatt6. 

—  The  power  returned  from  the  initial  wall  reflection  is  PrX  =  pPt  (or  PrX  dBm  =  PidBm+  Lp  dB ). 


6  There  is  a  “bookkeeping”  issue  here.  L  is  a  negative  quantity,  so  we  unintuitively  add  the  loss. 
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-  The  secondary  reflection  is  calculated  Pr2  =  p(rPt),  or  Pr2dBm  =  PidBm  +  LrdB  +  LpdB ,  in 
dB. 

—  Finally,  the  power  available  to  penetrate  the  wall  is  Pm  =  z1Pi  (or  Pm  dBm  =  PidBm+  2 LT  dB ). 


incident  signal 
P 


Pri  pPj  Pr2-pTPj 
signal  reflections 


Dielectric  medium  (wall) 


transmitted  signal 

p,= ^ 


>- 


Figure  38.  Schematic  of  reflection  and  transmission  coefficients. 


Numerical  example:  Consider  an  initial  pulse  power  Pt  of  1  watt  (30dBm),  with  sx—  1  and 
s2  =  10.  This  results  in  n:  =  1  and  n2  =  3.16.  Therefore,  p  =  0.269  and  t  =  (l  — /?)  =  0.73  . 
The  associated  dB  loses  are  Lp  =  -5.70  and  Lr  =  -1.37  . 


Table  5.  Reflection  losses— numerical  example. 


Quantity 

power  [watts] 

Power  [dBm] 

Initial  power,  Pi 

1  w 

30  dBm 

1st  reflection,  Pri 

pPi  =  0.269  x  Pi  =  0.269w 

Pi  +  LP  =  30dBm  -  5. 7dB  =  24.3dBm 

Power  in  medium,  Pm 

rPi  =  0.73  x  Pi  =  0.73w 

Pi  +  Lr  =  30dBm  -  1.37dB  =  28.6dBm 

2nd  reflection,  Pr2 

pPm  =  ppP< 

=  0.269  x  0.73  =  0.196w 

Pi  +  Lp  +  LT 

=  30dBm  -  5. 7dB  -  1.37dB  =  22.9dBm 

Transmitted  power,  Pt 

iPm  =  rPi  =  0.732  x  Pi  =  0.532w 

Pi  +  2Lr  =  30dBm  -2x  1.37dB  =  27.26dBm 

This  results  in  an  overall  reflection  loss  of  2.7dB  for  one-way  transmission  through  the  wall. 
Note  there  were  a  number  of  simplifications  to  obtain  this  value.  First,  we  considered  a  solid, 
homogeneous  dielectric  rather  than  the  cindercrete  wall  containing  air  pockets.  Second,  no 
wall  attenuation  was  accounted  for.  This  would  have  come  about  by  having  a  complex 
dielectric  £  —  s'  —  io  /  CO£0 ,  where  cris  the  conductivity  and  co  the  frequency.  The  dielectric 

constant,  £0  =  8.854x10 ~n  F  /  m  . 
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A  closer  examination  of  reflections  as  the  signal  reaches  a  discontinuity  provides  information 
about  the  structure  of  the  received  waveform.  For  example,  Figure  39  follows  the  path  of  a 
pulse  in  time  as  it  passes  through  a  wall,  reflecting  off  of  a  target. 

Rx/Tx  location 


Dielectric  medium  (wall) 


Figure  39.  Time-space  diagram  showing  signal  propagation  through  dielectric  material  (Not  to  scale). 

Moving  from  top  to  bottom,  the  time  delay  of  the  received  signal  increases.  Anything  fixed  in 
space  such  as  the  wall  or  a  stationary  target  traces  out  a  vertical  strip  as  shown.  The  lines  with 
arrows  show  the  various  paths  that  can  be  taken,  where  the  line  slope  indicates  velocity — EM 
waves  moving  into  the  dielectric  slow  down.  In  principle,  there  will  be  an  infinite  number  of 
reflections,  but  since  higher  orders  are  severely  attenuated,  we  only  consider  the  first  few 
returns.  The  reflection  and  transmission  losses  are  indicated  on  the  diagram  similar  to  that  of 
Figure  38. 
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A  more  accurate  representation  of  the  time  dependence  comes  from  properly  accounting  for 
propagation  distances  and  velocity.  This  was  carried  out  numerically,  and  the  signal  power  vs. 
time  delay  is  shown  in  Figure  40.  Note  the  order  in  which  pulses  were  received  differs  from 
what  is  indicated  in  Figure  39.  This  is  expected  since  that  diagram  was  only  used  to  track  the 

topology  of  paths  and  derive  terms  such  as,  r2 cr 0 per 0r2 Pt .  The  time  and  range  scales  in 
Figure  39  were  altered  simply  to  display  the  path  information  more  clearly. 


Figure  40.  Time  delay  vs  signal  power  after  multiple  reflections  from  a  dielectric  wall. 
Numbered  labels  refer  to  return  paths  as  indicated  in  Figure  39. 

This  is  an  useful  result  as  it  shows  the  relative  effects  of  reflection  losses  and  propagation 
delays  and  how  they  affect  the  overall  waveform  shape  and  attendant  pulse  shifting  due  to 
through  wall  propagation.  This  provides  a  relatively  easy  and  quick  way  to  investigate  the 
influence  of  dielectric  constant  and  the  target  RCS.  This  configuration  is  for  the  monostatic 
case,  but  could  be  modified  to  examine  bistatic  situations. 
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8.  Wall  Clutter 


As  discussed  in  the  previous  section,  building  materials  act  sources  for  multiple  reflections.  In 
through-wall  imaging  applications,  this  becomes  an  important  factor  for  image  quality.  Here, 
we  consider  a  wall,  with  a  characteristic  RCS,  <r0 .  As  the  wall  is  an  extended  structure  and 
not  a  point  scatterer,  it  is  convenient  to  express  the  RCS  per  unit  area  of  the  surface  as 
cr°  =  <r0  j A  .  The  symbol  cr°  is  referred  to  as  the  normalized  RCS,  or  backscatter  coefficient. 

The  magnitude  of  <j°  is  strongly  dependent  on  the  angle  of  incidence.  The  range  of  incident 
angles  expected  for  through-wall  application  is  relatively  small,  and  the  limiting  case  of  zero 
incidence  angle  will  be  the  dominant  source  of  clutter  returns.  If  the  wall  were  to  reflect 

isotropically,  then  cr0  would  be  unity.  With  isotropic  reflection  in  the  half-space,  o 0  becomes 
2  (3dB).  A  relatively  smooth  surface  (Rayleigh  criterion)  exhibits  a  strong  specular  reflection 
( <T°  >  2 ),  and  is  dominant  for  normal  incidence  [22]. 


8.1  Signal-to-Clutter  ratio 

The  following  treatment  is  adapted  from  some  of  the  literature  in  altimeter  radar  [23].  To 
examine  the  signal-to-clutter  ratio,  we  consider  the  received  power  from  a  target  of  RCS  <Jtgt , 

and  an  arbitrary  extended  structure  with  RCS  ac ,  where  the  subscript  denotes  “clutter”.  From 
the  standard  radar  equation,  we  write  the  signal  S,  and  clutter  C,  received  powers  as 


_  PrGAealgt 

Mr4 

_  PTGAec rc 
(4  x)3R4 


(26) 


(27) 


respectively,  where  Ae  is  the  effective  area  of  the  radar  antenna.  Substituting  cr°  =  cr c/  Ac  , 
the  signal-to-clutter  ratio  is 


S  =  = 

C  <7C  Aca° 


(28) 


The  value  Ac  is  the  area  of  the  radar  resolution  cell  and  depends  on  the  antenna  beamwidth 
and  range.  We  consider  nearly  perpendicular  incidence  angle,  and  the  instantaneous  field  of 
view,  or  illumination  area,  is  easily  obtained.  Generally,  an  antenna  is  characterized  by  its 
horizontal  and  vertical  beamwidths  0H ,  9V  ,  but  assuming  azimuthal  symmetry  yields 

9  H  =  9V  =  9zdB  .  At  range  7?,  such  an  antenna  will  illuminate  an  area  of  diameter  D  as 
detailed  in  the  following  diagram: 
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The  area  illuminated  is  therefore  Ac  =  7r(D / 2)2  =  7rR202dB  / 4 ,  and  the  resulting  signal  to 
clutter  (S/C)  ratio  is 


S_ 

C 


tgt 


Act0 


4  a 


tgt 


nR2eldBG° 


(29) 


Note  that  the  range  dependence  scales  as  HR2  rather  than  HR4  as  in  the  noise-limited  case. 
Also,  the  transmit  power  PT  does  not  appear  explicitly  although  it  remains  an  important 
parameter:  increased  transmit  power  increases  the  received  signal,  but  with  a  simultaneous 
increase  to  clutter  power,  yielding  no  overall  improvement. 


8.2  Transient  response 

We  next  wish  to  determine  the  power  returned  from  an  area  element  of  the  wall.  Rather  than 
simply  examining  a  point  scattering  process,  we  integrate  all  area  elements  providing  returns 
at  the  same  instant.  There  are  three  factors  that  affect  this-the  relative  geometry  between  radar 
and  scattering  surface,  the  antenna  pattern  on  the  wall  surface,  and  the  signal  duration. 

For  the  purpose  of  this  analysis,  the  geometry  remains  rather  simple.  The  source  of  clutter  is  a 
flat  surface  at  a  distance  R  from  the  antenna.  The  antenna  pattern  is  represented  by  an  angle- 
dependent  gain,  G{6\  Referring  to  the  experimental  pattern  (Figure  22),  a  plausible 

functional  form  is  G(o)  =  G(o)cosfl  6 ,  where  the  parameter  a  controls  antenna  beamwidth- 
the  larger  the  value,  the  narrower  the  beam.  The  pulse  duration  is  denoted  x. 

Further,  we  consider  the  backscatter  coefficient  to  be  of  the  form  cf°(#)  =  cr°(o)cos^  6.  The 
parameter  b  adjusts  surface  roughness  (for  a  smooth  surface,  b  is  large). 

The  received  power  becomes 


_  PTG2  (9)A2 <t° (o) 

Mr4 


(30) 
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We  next  determine  how  the  received  power  changes  in  time  given  that  both  the  antenna 
pattern  and  backscattering  coefficient  display  an  angular  dependence.  The  time  interval  is 
limited  to  small  multiples  of  the  pulse  width,  say  7  Ox.  For  an  infinitesimal  area  of  radius  z  and 
width  dz ,  at  time  t  we  have  (refer  to  Figure  41) 


Figure  41.  Geometry  used  to  calculate  clutter  pulse  returns. 

Even  though  we  consider  normal  incidence,  there  still  exists  a  narrow  range  of  angles  over 
which  the  pulse  interacts  with  the  wall.  The  first  task  therefore,  is  to  relate  the  incident  angle, 
6  to  the  radius  of  the  illuminated  area,  z.  In  the  following  diagram,  distance  perpendicular  to 
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the  wall  is  denoted  R0  (at  time  t  =  0 ).  At  a  later  time  t ,  the  area  of  illumination  has  increased, 
along  with  the  propagation  delay,  which  is  now  ct/2.  The  radius  of  the  illumination  disk  is  zt. 


Note  that  we  can  factor  out  the  cos2  9  term  using  (cos  Q^aJrb  -  (cos2  6 )  .  The  following 

relations  are  derived  from  the  above  diagram,  starting  with  cos  6  =  R0/R: 


cos  6  = 


Rt 


D  2  .  2 

R0  +zt 


1  +  zt  Rl 


R2=z2t+Rl 


z 


t 


CtRr 


1  + 


ct 

4 R, 


(32) 

(33) 

(34) 


These  are  required  because  the  integral  we  wish  to  solve  does  not  have  the  variable  z 
appearing  explicitly  in  the  integrand.  Substitution  into  the  following  equation,  carried  over 
from  the  previous  page, 


pA>) 


PTA2G-{0)a',(0)  !•<  (cos’  0) 

Jo 


{a+b/  2 


32;r 


R q 


-zdz 


(35) 


results  in 
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PM 


PTfG2(oy(o) «,  (l/(i  +  zV^))' 


M+b/  2 


12k2 

PtA2G2( 0)cr°(o)  (•-, 


-zJz 


12k2  Rl 


r 

Jo 


(l  +  z2/R2) 


■  dz 


(36) 


where  A:  =  2+a+Z>/2.  The  variable  z  now  appears  explicitly  in  the  integrand.  Performing  the 
integration,  received  power  from  the  wall  is 


pA‘)= 


PtA2G2  (o)cr°  (o)  1 


12k2  Rl 


2(k-\) 

Writing  out  the  explicit  time  dependence  of  zt,  we  have 


f  z2l 

r 

1- 

\ 

t 

R-  I 

1 

'v  iyo  J 

J 

z2  r 
1+iT  = 

K 


i+ 


ct 


\2 


v  2*0  j 


(37) 


(38) 


(]  *r22G2(Q)rT0(0)  1  f  1 

12k2  R2  2(A:-l)[  (l  +  ct/2R0  )2(M) 


(39) 


The  time  reference  t  =  0,  is  set  when  the  leading  edge  of  the  transmit  pulse  is  received  after 
being  reflected  from  the  nearest  part  of  the  wall.  The  area  of  the  wall  is  illuminated  for  t  <  r, 
at  which  point  the  center  of  the  illuminated  area  opens  up  forming  a  ring  of  illumination. 

The  return  signal  can  be  separated  into  two  regimes,  where  the  first  is  increasing  power  as  the 
direct  return  impinges  on  the  receiver.  The  duration  of  this  is  equal  to  the  pulse  width,  r.  After 
that,  the  signal  power  drops  off,  but  at  a  reduced  rate  depending  on  the  smoothness  of  the 
wall.  The  rougher  the  surface,  the  greater  will  be  the  range  of  scattering  angles.  This  translates 
into  a  slower  drop-off  for  the  received  power.  Beyond  time  r,  the  received  power  is  due  to  an 
expanding  ring  on  the  wall  surface.  To  obtain  this  contribution,  consider  the  following 

(0  =  Pr  {^)~Pr  (/  —  T)  (  40  ) 

for  t  >  r  .  At  t  =  r ,  we  have  PF  (t )  =  PR  (t )  -  PR  (o) . 

Recall  that 


zt 

PR  ( t )  =  2k  z  dz  ( 41 ) 

o 
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Where  we  are  integrating  over  the  shaded  area  in  the  following  diagram: 


zt 

PF(t)=  |  PR27izdz 

Zt-T 

0  Zt 

=  J  Pr2k z dz  +  ^ Pr2k z dz  (42) 

zt-r  0 

Zt  Zt-T 

=  ^PR27rzdz-  | PR27izdz 
0  0 


This  last  expression  is  rewritten  as 

(0  ~  Pr  (0  Pr  ^")  *  (  43  ) 

An  example  of  the  pulse  tail  duration  is  given  in  Figure  42,  where  we  have  set  R0  =  10m,  r  =  1 
ns,  and  the  antenna  beamwidth  factor,  a  =  2. 


Figure  42.  Time  signal  of  pulse  returns. 
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Having  worked  through  the  analytic  details,  it  is  relatively  straight-forward  to  implement  this 
in  a  few  lines  of  code. 

The  plots  in  Figure  42  were  produced  using  the  Matlab  routine  in  Table  6. 

Table  6.  Matlab  code  for  wall  clutter  return  pulse. 


% - Wall  Clutter - 

a  =  2;  %  antenna  beamwidth  factor 

b  =  10;  %  surface  roughness  factor 

%  (smooth  surface;  b  large) 
k  =  2  +  a  +  b/2; 
gamma  =  2*(k-l); 
tau  =  1 .0;  %  pulse  width  [ns] 

c  -  0.2998;  %  [m/ns] 

d  =  10;  %  range  from  wall  [m] 

tmax  =10;  %  maximum  delay  [ns] 

Xi  =  c/2/d;  %  [1/ns] 

A0  =  1 ;  %  amplitude 


del_t  =  0.01; 

tl  =  0:del_t:tau; 

t2  =  tau+del_t:del_t:tmax; 


%  P_rise:  power  received  as  direct 
%  pulse  returns  from  wall  extends 
%  to  t  =  tau,  pulse  duration 
P_rise  =  A0*(1  -  (1  +  Xi*tl).A-gamma); 

%  P_fall:  power  fall-off  due  to 
%  finite  range  of  incident  angles 
P  fall  =  A0*(  (1  -  (1  +  Xi*t2).A-gamma). . . 

-  (1  -  (1  +  Xi*(t2-tau)).A-gamma) ); 

%  normalized  plots 
plot(t  1  ,P_rise/max(P_rise)) 
hold  on,plot(t2,P_fall/max(P_rise); 
ylabel(’P_R/(P_R)_  {max} f) 
xlabel(?t  [ns]?) 
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9.  Application to avalanche detection 


One  particular  application  of  UWB  impulse  radar  is  to  assist  in  search  and  rescue  operations 
for  avalanche  victims.  An  obvious  advantage  of  this  technology  is  the  ability  to  generate  high- 
accuracy  range  data,  while  at  the  same  time  achieving  good  EM  penetration  with  center 
frequencies  less  than  two  Gigahertz. 

As  mentioned  in  the  introduction,  portions  of  this  work  were  sponsored  by  the  Canadian 
National  Search  &  Rescue  Secretariat  (NSS)  to  investigate  the  feasibility  of  using  short-pulse 
UWB  radar  to  detect  objects  buried  in  snow.  This  is  a  Canadian  multi-department  venture 
including  Parks  Canada,  the  NSS,  and  the  Department  of  National  Defence. 

The  experimental  configuration  [24]  consisted  of  a  large  mound  of  snow,  approximately  32 
m2,  piled  up  next  to  a  small  building  that  housed  the  equipment.  The  snow-building  boundary 
was  composed  of  wood  (plywood  and  2  x  4?s). 

9.1  Range  processing 

In  any  two-  or  three-dimensional  radar  imaging  system,  the  basic  element  of  target 
information  is  a  time-delayed  impulse  response.  It  is  therefore  convenient  to  first  consider 
one-dimensional  range  data  to  assist  with  feasibility  studies  for  such  a  system. 

The  objective  of  our  experimental  work  was  to  measure  radar  returns  from  various  objects 
buried  in  a  mound  of  snow,  examining  the  approximate  SNR  of  the  target  response. 

The  results  shown  here  are  for  a  single  transmitted  pulse.  No  signal  enhancement  was 
obtained  via  pulse  integration.  Data  processing  consisted  of  application  of  a  bandpass  filter 
followed  by  background  subtraction.  Here,  background  subtraction  reduced  dynamic  range  — 
the  initial  impulse  from  the  transmit  antenna  was  removed,  as  well  as  reflections  from  a  tunnel 
in  the  snow  (used  to  place  targets).  This  environment,  in  the  absence  of  a  target,  formed  the 
reference  background. 

Figure  43  shows  returns  from  the  snowbank,  where  various  objects  were  placed.  The  wall, 
tunnel  and  snowbank  acted  as  background  reference,  and  was  subtracted  from  each 
subsequent  data  set.  The  notation  “AVH”  denotes  the  type  of  pulse  generator  used  in  the 
experiment  (Avtech  Pulser  model  #  AVH-HV1-C-T1).  Range  values  need  to  be  calibrated, 
but  the  targets  appear  at  a  relative  range  of  ~2.3m,  which  is  the  approximate  radial  distance  to 
targets. 
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Range  [m] 

Figure  43.  Measured  data  after  background  subtraction  for  different  buried  targets. 
This  data  set  consists  of: 

a.  Reference  return  with  snowbank  and  tunnel 

(All  subsequent  plots  have  this  reference  background  subtracted). 

b.  Radar  return  from  human  phantom  filled  with  saline 

c.  Radar  return  from  a  metallic  sphere. 

d.  Radar  return  from  metallic  cylinder. 

e.  Radar  return  from  a  wooden  block. 


Target  returns  from  the  metallic  objects  display  the  largest  values,  while  the  wooden  block 
produced  the  least  reflections.  Returns  due  to  the  human  phantom  lie  between  these  two 
extremes  and  is  clearly  visible. 


9.2  SAR  processing 

The  system  configuration  (Figure  44)  consists  of  one  transmit  antenna  at  coordinates  (x=0.0, 
y=0.58,  z= 0.64)m.  The  receive  antenna  is  moved  along  the  x-axis,  forming  a  synthetic  array. 
In  the  experiment,  there  were  five  elements  in  the  array  at  a  spacing  of  0.5m,  generating  an 
effective  length  of  2.0m.  They-  andz-  coordinates  were  fixed  at  zero.  The  imaging  area  has 
dimensions  of  4m  (azimuth)  X  5m  (range).  The  ’’wall”  in  this  experiment  is  composed  of 
plywood  sheets  acting  as  a  boundary  between  the  lab  and  the  snowbank.  The  plywood  is 
located  0.8m  in  front  of  the  Rx  array,  and  0.2m  in  front  of  the  Tx  element. 
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Figure  44.  Schematic  showing  the  experimental  configuration. 

9.2.1  Simulations 

Sets  of  simulations  corresponding  to  the  experimental  configuration  were  carried  out,  to 
determine  if  the  qualitative  features  of  SAR  images  could  be  reproduced.  Within  these 
simulations,  two  point  reflectors  were  located  at  (x=-0.5,  y= 2.5,  z=0.0)m;  <j=  1.0  and  (x=0.0, 
y= 4.0,  z=0.0)m;  cr=  2.0.  Here,  the  parameter  o denotes  the  reflection  coefficient. 

The  first  plot  (Figure  45a)  shows  the  image  reconstruction  using  the  five  element  positions  as 
described  previously,  only  now  for  a  monostatic  configuration.  This  means  that  Rx/Tx 
antennas  are  co-located  and  move  along  the  synthetic  aperture  as  one  unit.  Because  of  the 
greater  angular  diversity  inherent  in  such  a  system,  the  resolving  capability  is  better  than  that 
of  the  bistatic  case  (Figure  45b).  For  example,  resolution  can  be  characterized  in  terms  of  full- 
width  at  half-maximum  (FWHM),  and  for  the  two  cases,  we  obtain  the  following  ( In  both 
cases,  range  =  2.5m,  synthetic  aperture  length  =  2m ): 

•  Monostatic:  FWHM  =  1.1  m; 

Bistatic:  FWHM  =1.7  m. 

This  yields  a  ratio  of  0.65,  close  to  the  theoretical  value  of  one-half  obtained  from  comparison 
of  a  synthetic  aperture  vs.  a  physical  array. 
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Figure  45.  Simulated  data,  a)  Monostatic  configuration,  b)  Bistatic  configuration. 


Figure  46  shows  a  case  where  a  “wall”  has  been  simulated  by  generating  a  line  of  point 
reflectors.  This  is  again  for  the  bistatic  case.  The  returns  from  desired  targets  are  clearly 
present,  but  are  severely  reduced  in  relative  magnitude  compared  with  the  wall  returns.  In  the 
absence  of  the  wall,  this  would  revert  back  to  the  image  depicted  in  Figure  45b.  Somewhat 
artificially,  the  wall  reflection  coefficient  has  been  set  arbitrarily  low  at  <j=  0.5,  but  even  this 
shows  how  drastic  the  effect  can  be.  This  does  not  account  for  multipath,  nor  does  it  mimic 
propagation  through  a  dielectric  layer,  but  it  does  capture  the  effects  of  increased  dynamic 
range. 


Figure  46.  Increased  dynamic  range  accompanying  the  presence  of  a  “wall”  of  point  reflectors. 
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9.2.2  Experimental  data 


The  collected  data  was  processed  using  the  standard  method  of  time-domain  back  projection 
to  focus  returns  based  on  time  delay.  The  kernel  of  this  algorithm  is  the  same  as  that  used  in 
simulations.  Immediately  we  can  see  two  of  the  same  characteristics  highlighted  in  simulated 
images:  dynamic  range  and  smearing  out  of  target  responses  due  to  the  bistatic  configuration. 


Figure  47  shows  the  processed  data  with  wall  returns.  As  expected,  this  dominates  the  image 
and  the  desired  target  (human  phantom)  is  not  at  all  obviously  visible. 


Range  [m] 


Figure  47.  No  background  subtraction.  Wall  returns  dominate  the  image. 


To  better  view  the  target,  the  same  scene  without  the  human  phantom  was  recorded  and 
subsequently  subtracted  from  the  original  image.  The  result  is  shown  in  Figure  48. 


Figure  48.  With  background  subtraction.  Wall  returns  negated,  two  strong  reflectors  beyond  human 

phantom  dominate  the  image. 
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10.  Discussion  and  Conclusions 


UWB  impulse  radar  is  inherently  a  noise-limited  technology.  While  sub-nanosecond  pulses 
achieve  good  range  resolution,  the  resultant  extreme  bandwidth  extracts  a  toll  on  system 
performance  in  terms  of  noise  power  entering  the  receiver.  As  is  typical  of  impulse  radar 
returns  that  lack  the  filtering  commonplace  in  narrowband  systems,  individual  returns  are 
severely  impaired  by  thermal  noise,  clutter  and  external  interference.  Data  must  be  processed 
efficiently  to  yield  an  acceptable  signal  to  noise  ratio  (SNR),  and  to  enhance  structural  details 
of  the  return  pulse.  This  is  achieved  by  a  combination  of  ensemble  averaging  to  reduce 
thermal  fluctuations  and  the  newly  developed  Background  Noise  Conditioning  (BNC) 
method,  a  statistical  signal  processing  technique  used  to  remove  isolated  spectral  interferers. 
BNC  was  developed  at  DRDC  Ottawa  to  periodically  examine  spectral  noise  properties,  and 
automatically  configure  an  appropriate  notch  filter. 

For  the  radar  data  presented,  averaging  over  a  relatively  small  number  of  pulses  is  insufficient 
to  raise  the  desired  returns  appreciably  above  the  noise  floor.  Typical  SNR  values  at  N  =  1  are 
around  18  to  19  dB.  Beyond  ten  pulses,  SNR  improves  approximately  linearly.  Integration  of 
about  80  pulses  results  in  a  SNR  maximum  of  about  35  dB ,  with  an  associated  improvement 
of  -15  dB  over  the  untreated  data  (that  is,  for  situations  where  BNC  is  not  applied).  For  values 
N  >  80 ,  the  averaging  takes  place  over  a  long  enough  time  period  that  returns  become  blurred 
due  to  target  motion.  This  effectively  defines  the  upper  limit  of  N  that  can  be  achieved 
without  additional  processing  to  focus  the  returns. 

Using  the  space-time  backprojection  technique,  reconstructed  data  clearly  shows  the 
oscillatory  motion  of  the  target,  but  the  signal  vanishes  as  the  radial  velocity  approaches  zero 
at  the  two  extremes  of  motion.  Again  using  experimental  data,  we  were  able  to  verify  the 
minimum  detectable  velocity  for  the  system,  depending  on  sampling  increment,  PRF  and  the 
number  of  pulses  integrated,  N.  This  information  can  be  employed  to  construct  a  non-doppler 
MTI  filter  for  target  velocity  discrimination.  In  turn,  this  could  be  used  to  further  enhance 
SNR  by  time-aligning  the  integrated  return  pulses. 

SNR  is  the  metric  that  is  most  readily  identifiable  in  characterizing  radar  detection.  This, 
together  with  the  radar  range  equation,  becomes  the  basis  for  radar  performance  prediction. 
This  model  was  employed  semi-empirically,  using  parameter  values  that  could  either  be 
measured  directly  or  elucidated  from  available  UWB  experimental  data,  such  as  antenna  gain, 
noise  figure  and  background  noise  power.  Independently,  SNR  and  the  path  loss  exponent 
were  obtained  directly  from  recorded  data.  Agreement  with  theoretical  trends  was  quite  good, 
although  a  recurrent  issue  is  that  close  range  data  appears  to  have  lower  SNR  than  predicted. 

A  tentative  conclusion  is  that  at  close  ranges,  the  system  is  unable  to  properly  gate  out  the 
clutter  returns,  artificially  raising  the  noise  portion  of  the  signal  to  noise  ratio.  This  behaviour 
is  also  observed  in  the  reconstructed  ISAR  images. 

Part  of  this  report  is  devoted  to  understanding  and  characterizing  the  impulse  shape  and  its 
associated  spectrum.  Contrary  to  narrowband  systems,  bandwidth  is  determined  less  from 
pulse  width  than  rise  time.  For  example,  a  0.5ns  pulse  width  will  not  always  translate  to  a 
2GHz  bandwidth,  as  is  often  reported  in  the  literature.  For  the  experimental  system  employed 
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here,  the  instrumentation  produces  a  Gaussian  pulse,  where  the  pulse  width  and  rise  time  are 
not  independent  parameters  —  specifying  one  fixes  the  other.  Another  issue  related  to  proper 
characterization  of  a  system  is  how  the  pulse  width  is  defined.  Due  to  its  easy  experimental 
access,  we  chose  to  use  the  full  width  at  half  maximum  (FWHM)  as  the  metric.  From  this,  we 
had  a  center  frequency  of  380MHz,  3dB  bandwidth  of  220MHz,  and  an  effective  pulse  width 
of  1.3  ns  for  an  initial,  transmit  pulse  width  of  0.42ns. 

A  wall  or  similar  dielectric  boundary  plays  a  role  in  the  application  of  UWB  systems.  The  two 
major  concerns  are  attenuation  and  clutter.  Taking  a  series  of  measurements,  we  were  able  to 
determine  the  cindercrete  wall  induced  about  8  dB  attenuation,  while  the  same  wall 
thouroughly  wet  caused  only  an  additional  2dB  loss  (both  of  these  are  one-way 
measurements). 

In  through-wall  imaging,  the  building  material  acts  as  a  reflector,  with  its  own  characteristic 
RCS,  a  .  As  the  wall  is  an  extended  structure  and  not  a  point  scatterer,  it  is  convenient  to 
express  the  RCS  per  unit  area  referred  to  as  the  normalized  RCS,  or  backscatter  coefficient. 
Using  this  as  the  starting  point,  we  analyzed  how  the  pulse  would  reflect  off  the  wall,  as  a 
general  function  of  the  backscatter  coefficient.  An  important  consideration  is  that  even  though 
we  consider  normal  incidence,  there  still  exists  a  narrow  range  of  angles  over  which  the  pulse 
interacts  with  the  wall,  leading  to  extended  tails  on  the  reflected  pulses.  This  will  lead  to  a 
better  understanding  of  how  to  set  the  PRI  to  reduce  this  wall  clutter. 

Finally,  we  began  to  investigate  imaging  of  objects  buried  in  snow,  using  both  ID  ranging  and 
2D  SAR  techniques.  The  results  show  that  objects  can  be  located,  although  the  current 
resolution  is  such  that  it  would  make  identification  challenging.  A  major  concern  would  be 
processing  technique  itself,  where  the  ensemble  averaging  performs  pulse-to-pulse  subtraction 
to  minimize  thermal  noise.  This  is  a  form  of  area  MTI,  where  returns  from  stationary  targets 
are  partially  removed  in  the  processing  stage.  These  results  will  be  used  to  develop  algorithms 
for  extracting  stationary  targets  without  apriori  knowledge  of  the  background.  One  approach 
could  be  to  establish  a  data  bank  of  target  signatures  and  use  these  to  process  data  returns  via 
correlation  techniques. 
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